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SUMMARY 


The  methods  of  Hamilton  and  Lagrange  in  the  calculus  of 
variations  are  employed  to  optimize  planetocentric  orbital  transfers 
using  fixed-thrust  engines,  where  the  thrust  acceleration  <(<  local 
gravitational  acceleration.   The  equations  were  derived  for  minimum 
time  transfers.   However,  since  fixed-thrust  engines  are  used,  the 
minimum  fuel  problem  is  solved  simultaneously.   Numerical  solutions  are 
obtained  by  integrating  the  optimized  differential  equations  of  motion 
for  the  problem,  the  actual  integrations  being  performed  by  an  IBM  360/30 
computer.   The  numerical  integrations  result  in  optimal  spiral 
trajectories  representing  orbital  transfers  between  specified  coplanar, 
co-axial  ellipses,  with  the  optimal  thrust  angle  being  calculated  at  all 
points  along  the  trajectories.   For  an  orbital  transfer  at  the  planet 
Mars,  the  low  fixed-thrust  system  yields  a  significant  increase  in  the 
amount  of  payload  that  can  be  delivered  into  the  desired  final  orbit 
compared  with  a  transfer  accomplished  by  an  impulsive  (high  thrust)  chemical 
rocket  propulsion  system.   However,  the  time  required  to  achieve  the 
transfer  is  necessarily  longer  in  the  low-thrust  case.   It  is  then  necessary 
to  examine  the  trade-off  between  the  increase  in  payload  and  the  extra 
transfer  time  to  determine  whether  or  not  the  fixed-thrust  system  is 
desirable. 
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I.   INTRODUCTION 

As  we  explore  the  space  around  us,  we  look  for  ways  of  doing  things 
more  efficiently,  ways  to  save  money,  and  ways  to  put  heavier  payloads  into 
orbit  at  the  same  or  less  cost  than  is  presently  possible.   One  feasible  way 
of  accomplishing  this  is  through  the  use  of  solar-electro-static  or  nuclear- 
electro-static  propulsion  systems.   They  require  much  less  fuel  (weight- 
wise)  than  chemical  rocket  engines,  have  much  higher  specific  impulses,  and 
can  result  in  considerable  increases  in  final  payload  delivered.   They  do, 
however,  require  more  time  to  achieve  the  same  orbital  maneuvers  because  of 
the  much  lower  thrust  inherent  in  these  type  of  systems. 

Considerable  research  on  optimal  low  thrust  maneuvers  has  been  done 
in  recent  years.   A  large  part  of  this  work  has  been  done  by  T.  N.  Edelbaum. 
He  has  published  papers  which  contain  solutions  for  optimum  power-limited 
transfer  and  rendezvous  between  arbitrary  coplanar  and  coaxial  elliptic  and 
circular  orbits.   The  assumptions  involved  in  the  development  are  1)   the 
thrust  acceleration  is  much  smaller  than  the  local  gravitational  acceleration. 
2)   the  propulsion  system  is  power-limited,  i.e.,  it  operates  at  constant 
exhaust  power  with  a  variable  thrust  magnitude  inversely  proportional  to  the 
exhaust  velocity.   3)   the  direction  and  magnitude  of  the  thrust,  both  of 
which  are  assumed  to  be  completely  variable,  is  to  be  determined  as  a  function 
of  time  so  as  to  minimize  the  fuel  consumption  for  a  fixed  total  transfer  time, 
Edelbaum  has  also  developed  an  analytic  solution  for  fixed  thrust  systems,  but 
this  is  only  applicable  to  transfers  between  almost  circular  orbits,  for  small 
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changes  in  a   and  e  ,  and  for  large  changes  in  a   and   i   (orbital  plane 
inclination) (Reference  1,  2,  3). 

This  thesis  addresses  the  problem  of  optimizing  fixed,  low- 
thrust  orbital  transfers  with  large  changes  in  a   and  e  .   Based  upon  a 
Lagrange-Hamilton  formulation  developed  by  Dr.  M.  Handelsman  (Reference  4) 
the  numerical  solution  of  optimized  transfer  trajectories  is  performed  by 
computer.   The  solutions  involve  the  optimization  of  planetocentric  transfers 
between  coplanar,  coaxial  ellipses  or  circles  of  specified  characteristics, 
with  no  time  constraint.   Because  of  the  relatively  small  fuel  expenditure 
in  these  transfers,  it  is  possible  to  replace   "fixed  thrust"  with  "fixed 
thrust  acceleration"  throughout  the  entire  transfer  with  little  error.   Solar 
electric  engines  lose  a  certain  amount  of  thrust  as  their  distance  from  the 
sun  increases;  this  could  hardly  be  considered  constant  thrust  in  heliocentric 
space.   But,  since  the  problem  has  been  limited  to  a  planetocentric  study, 
the  thrust  of  the  solar  electric  engine  can  be  assumed  to  be  constant  within 
the  planetocentric  space. 

By  fixed  thrust,  it  is  meant  that  the  thrust  magnitude  is  constant 
and  that  the  engine  remains  turned  on  at  all  times  during  the  transfer.   The 
problem  then  remains,  where  to  point  the  thrust  at  all  times.   Now,  there  are 
many  factors  which  could  be  optimized  in  an  orbital  transfer  of  this  type. 
In  this  thesis  two  separate  phenomena  were  optimized.   First,  orbital  trans- 
fers were  optimized  for  minimum  time  using  fixed  thrust  engines.   Since  the 
thrust  is  fixed,  and,  therefore,  continuous,  then  a  minimum  time  transfer 
would  also  be  a  minimum  fuel  transfer,  the  amount  of  fuel  expended  being 
directly  proportional  to  the  length  of  time  the  engine  is  turned  on.   Secondly, 
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orbital  transfers  were  optimized  for  minimum   E   (eccentric  anomaly);  in 
other  words,  transfers  were  achieved  to  the  desired  orbit  in  the  least  amount 
of  revolutions.   The  significance  of  this  type  of  transfer  trajectory  is  that 
it  lends  considerable  insight  into  the  minimum  time,  minimum  fuel  problem. 

The  problem  formulation   involves  deriving  the  perturbation  equations 
of  motion  for  fixed  thrust  engines  of  low,  constant  magnitude,  and  then  opti- 
mizing them  for  minimum  time  and  minimum  E   using   Hamilton  ,  Lagrange,  and 
Mayer  techniques  in  the  calculus  of  variations.   (Reference  4).   After  being 
optimized  for  minimum  time  and  minimum  E  respectively,  the  differential 
equations  are  set  up  to  solve  for  the  optimum  thrust  angle  necessary  at  all 
times  to  achieve  the  optimized  transfers.   These  sets  of  equations  are  then  solved 
by  numerical  integration, in  this  case  on  an  IBM  360-30  computer  system.   An 
analytic  solution  is  not  available  at  present  because  of  the  inability  to  solve 
for  Langrangian  multipliers.   This  is  the  reason  for  the  numerical  approach. 
The  basic  assumptions  follows: 

1)  The  thrust  acceleration  is  much  less  than  the  local  gravitational  acceler- 
ation.  (A«g). 

2)  Continuous,  fixed-magnitude  thrust. 

3)  Relatively  small  total  fuel  consumption;  therefore,  for  the  sake  of  simpli- 
city, the  thrust  acceleration  is  assumed  to  be  essentially  constant  over  the 
transfer.   That  is,  constant  thrust  magnitude  is  taken  as  synonymous  with 
constant  thrust  acceleration.   (A  program  modification  to  incorporate  the 
changes  in  thrust  acceleration  due  to  fuel  depletion  is  simple  to  implement). 

A  description  of  the  problem  and  its  formulation  follows. 
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II.   DESCRIPTION  OF  PROBLEM 

An  initial  elliptic  orbit  is  assumed  about  a  planet  or  any  other 
body  which  has  an  inverse-square  gravitational  field.   The  ellipse  has 
specifications   (aQ,  e  ),  where   a    is  the  initial  semi-major  axis  and 
e    is  the  initial  eccentricity.   The  problem  then  is  to  transfer  optimally 
to  a  desired  final  coaxial,  coplanar  ellipse  with  specifications   (a  ,  e  ). 
See  figure  1.   This  is  accomplished  by  turning  on  the  engine  at  the  initial 
ellipse  periapse  and  letting  the  optimally-directed  fixed  thrust  program  guide 
the  vehicle  along  a  spiral  trajectory  until  the  final  ellipse  is  obtained. 
The  engine  is  then  shut  off  and  the  spacecraft  is  in  its  new  orbit,  the  trans- 
fer having  been  conducted  in  minimum  time  and  fuel  expenditure  or  in  the  least 
amount  of  revolutions,  depending  upon  which  program  is  used.   The  numerical 
results  show  that  the  minimum  time  transfers  always  take  more  revolutions 
than  the  minimum  E   transfers,  and  the  minimum  E   transfers  always  take  more 
time  than  the  minimum  time  transfers.   This  is  as  expected  and  lends  consider- 
able support  to  the  calculations. 

Each  particular  transfer  is  determined  by  the  selection  of  the  initial 
Lagrangian  multiplers.   Given  an  initial  ellipse,  different  sets  of  inital 
Lagrangian  multipliers  will  cause  the  vehicle  to  transfer  to  different  ellipses. 
Since  we  have  no  analytic  way  of  picking  these  multipliers  at  present,  many 
cases  are  run  to  cover  a  wide  range  of  transfers  from  given  initial  orbits. 

The  formulation  of  both  the  minimum  time  and  minimum  E   equations 
follows  as  the  first  phase  of  the  problem  solution. 
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III.   LAGRANGE  FORMULATION  IN  VARIABLE   E 
A.   Orbital  Perturbation  Equations 

The  general  orbital  perturbation  equations  in  terms  of  radial 
tangential  thrust  forces,   S   and   T  respectively,  shown  in  Figure  2, 
are  as  follows:  (Reference  5) 

Some  useful  relations  follow: 


S 


m  (VEHICLE) 


FIGURE  2. 
Vehicle  Force  Vectors 


'■■  &  ■  -hr*  (sesin£  + 

n\l-e 
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(1) 


e   = 


de 

dt 


\KZ2 


n  a 


S  s  in  f  +  T(cos   E  +  cos    f ) 


(2) 


Using: 


:os   E     +      cos    f      = 


JL  (  J.  -  i) 

e     v    r  a  ' 


p  =  a(l-e  ) 


f  =   true  anomaly   =   0 


E  =  eccentric  anomaly 


Equations  (1)  and  (2)  become: 


a   = 


2ejsin   f}s               2a  Vl-e2 
i =s  +  


n  Vl-e 


r  n 


(3) 


yJT2      ■ 

Vl-e        sin 


na 


*-     .+^ 


2 

a   ner 


a    (1-e    )    -    r 


(4) 


2   3 

na        =y=M=k      =     GM 


(5) 


n   = 


Vk/a3        =      2   tt      / 


(period) 


(6) 


S    =      FNORMAL/m      =      F    sina     / 


(7) 


T=      FTANG./m      =      FCOSQ/ 


(8) 


For   an  ellipse 


r   =     *    (1-e    ) 
1+e   cos    f 


(9) 


Through  substitution  and  simplification  equation  (3)  yields: 


2      I 


VT7 


:# 


e  sin  8  +  sin«  +  (1  +  e  cos  0)  cosa 


(10) 
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Similarly,    equation    (4)    yields 


k  -    n/^72   ^yfl 


m     *  k 


_       .  2    cos    9  +   e(l   +   cos      8)    cos  a 

■  in  0   sin  a    +  : — * — « 

1  -f   e    cos    0 


(ID 


Next,    it    is   desired    to  write   equations    (10)    and    (11)    in   terms 
of   the   eccentric   anomaly     E.      The    following  relations   between     E      and      8 
are   used    (Reference   6): 


1  +  e   cos   0 
sin  0      =     — t — -  9 sin     E 

\/l-e 


(12) 


1  +  e   cos   0     = 


1-e 


1-e   cos   E 


(13) 


.     .      sin   0   = 


\l-e      sin  I 
1-e   cos   E 


(14) 


cos    0    ■ 


cos   E    -    e 
1-e   cos   E 


(15) 


r   = 


1+e    cos    0 


ad-e    ) 

"      1^ 

1-e   cos   E 


a(l-e   cos   E) 


(16) 


1+e   cos   0 


1-e 

1-e   cos   E 


(17) 


sin  E   =  sin  0 


d0   = 


V^? 


1-e    cos   E 


dE 


(18) 


(19) 
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D.   Comparison  with  Edelbaum's  Equations 

Substituting  equations  (12-19),  equations  (10)  and  (11)  can  be 
written  in  terms  of  E: 


2F 


a     (1-e  cos  E) 


e  sina  sin  E  +  cos 


oVl-e 


(20) 


e   = 


(1-e  cos  E)m 


2  I       2  2 

(1-e  )  sin  E  sina  +  \l-e 


(2  cos  E  -  e  cos   E-e)  cos  a 

(21)- 


Equations  (20)  and  (21)  agree  with  Edelbaum's  equations  for 
low  thrust  perturbation  (Reference  1). 
Since 


dt   =  dE 


Vv  (i- 


e  cos  E) 


(22) 


equations  (20)  and  (21)  become 


da     2Aa' 
dE       k 


[ 


e  sin  E  sin 


a   +  \l-e    COS  a 


(23) 


de     Aa_ 
dE      k 


2  I       2  2 

(1-e  )  sin  E  sina  +  Vl-e    (2  cos  E-e  cos  E-e)  cos  a 


(24) 


Since  we  are  optimizing  the  thrust  angle  for  minimum  flight  time 
the  following  integral  is  minimized: 

E, 


(1-e  cos  E)  dE 


(25) 
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Th  is    is    of    the    form    (Reference   4,    7) 


y«0  av  x2. 


E 

t£   = 


where: 


fo   =    /?     =\a  /k    (1-e   cos  E)  (27) 

Xx  =  a;   X2  -  e  (28) 

C.   Employing  the  Hamiltonian 

To  find  the  minimum  extremal,  we  write  the  Hamiltonian  for  this  system  of 
equations , 

2 
H  -  -  fQ  +  £    X.  (E)  f.   ,  (29) 

i=l 


where 


dxi       H, 
f   =  — i-  =  -^  (30) 

rl    dE       dE  l   ' 


2    de 
f2  "  dE   "  dE       '  (31) 


The  independent  variable  is  the  eccentric  anomaly  E.   Substituting 
into  H: 

H  -  "  V^k  ^  COS  E)  +   Xa  ol  +  Xe  I  (32) 
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^  ,  da       de 

Substituting  for  —  and  — 

ah       aE 


H  =   "V^/k   (1"e  C°S  E)  +  Xa"^ 


e  sin  E  sin  a    + 


V 


1-e   cos  a 


+  A 


Aa 


e   k 


(1-e  )  sin  E  sin  a  + 


V 


2  2 

1-e   (2  cos  E  -  e  cos   E  -  e)  cos  a 


(33) 


D.   Euler-Lagrange  Equations  for  Adjoint  Variables 


The  following  are  the  Euler-Lagrange  equations  for  the  adjoint  variables  a 


dA. 

l 

dE 


dn 
dx. 


(i  =  1,2) 


(34) 


Expanding  equation  (34) 


dA 


dE 


(A 


df- 


d*< 


+  A 


)  +    T  <  ^fi   +    A    f9) 


30 


a    d  a      '    "e       da      '  $     v  "a'l         Ae   2'       da 


(35) 


dX 


dE 


afi  3f2  l  dp 

<     ^"71     +    Xe     ~  )+7    <Aafl  +    A  eV  ~ 


(36) 


E.   Transversality  Condition 


Next,  the  transversality  condition, 


HdE  +  A  .  dX. 
l    l 


=  0 


(37) 
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which  is  to  be  evaluated  using  the  boundary  conditions 


t=0:    a=a,   e=e,   E   =0  (36) 

o        o    o  v  ' 


t  =  t  :   a  =  a  ,   e  =  e  ,   Ef  =  OPEN  (39) 

The  choice   E   =  0   is  made  for  convenience;  this  is  the  periapse  point,  the 
usual  choice  for  the  initial  retro-impulse  into  a  bound  planetocentric  elliptic 
orbit  from  an  incoming  hyperbolic  trajectory. 

Combining  equations  (37),  (38),  and  (39), 

Aafl  +  Aef2  "  fi     =  *  =  °  (at  t=tf}  (40) 

F.   Control  Equation  for  Optimal  Thrust  Angle 

Equation  (40)  is  a  necessary  condition  for  a  transfer  to  be  optimized  for 

minimum  time  (and  fuel).   For  the  optimal  thrust  angle  a    , 


3h  3fl         df2 

±H  =  0   =   A   — -  +  A   T^  (41) 

d<*  a  da  e      d  a 


which  becomes 

A    2Aa3 
a 


0      = [e    sin  E   cos  a       -        ^1-e      sin   a   ]      + 


^l-e      sin  a 


.    2 
A    Aa 
e 


2  I       21  2 

f(l-e    )    sin  E   cos  a    -     Jl-e      (2   cos   E   -   e   cos      E   -   e)    s in  a]  (42) 
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Divide  through  by  cos  «  and  solve  for  tan  a    : 

2 
2a  A   e  sin  E  +  A   (1  -e  )  sin  E  .,„. 
a e        (43) 

*   a    -yl-e        +    A  e  \l-e  (2   cos  E    -   e   cos      E    -   e) 

To   simplify   the   equation   for   computer   programming   purposes    let: 

X     =     e   sin  E  (44) 

Z     =      (1-e2)    sin  E  (45) 


=-    >P 


e"  (46) 


-      J7" 
W     -      >(1- 


l 
e*"    (2   cos   E   -   e   cos*"  E   -   e)  (47) 


N     z     2a    A       X  +    A       Z  (48) 

a  e 


D     z     2a   A        Y  +   A       W  (49) 

a  e 


Substituting  into  equation  (41) 


2a  *  a  X  +  A  e  Z  _  J. 
tan  a  =  ~2^~X       Y  +  A   W  =   D 


(50) 


This  is  the  control  equation  for  the  optimal  thrust  angle  a   .      The  step  variable 
of  integration  is   E.   Time  of  transfer   (t)   is  determined  by  integration  of 
equation  (22)  which  can  be  written: 


dt  _  JL_ 

di  "  p  (51) 
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Thus,  we  now  have  sufficient  equations  for  numerical  solutions  for 

different  transfers,  depending  on  the  values  assigned  to  \  and    \ 

o  o 

The  problem  is  next  solved  using  the  Mayer  formulation  to  develop  an 

alternate  but  equivalent  set  of  equations. 
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IV.   MAYER  FORMULATION  IN  VARIABLE   t 

The  Mayer  formulation  of  the  problem  proceeds  as  follows 


dX. 

-   f,(X.,  u    t)  (52) 


dt       iv"i'    j 


Where: 


Xx  =  a,  X2  =  e,  X3  =  E,  ux  =  a  (53) 


The  independent  variable  now  is   t   rather  than  E   as  in  the 
Lagrange  formulation;  here  we  minimize  the  value  of   t  ,  which  will  be 
called   G.   Thus: 


G  -  t£  (54) 


A.   Hamiltonian  and  Euler-Lagrange  Equations 

The  Hamiltonian  is  different  and  does  not  contain   f 


3 
H  =    £     A  .(t)  f.  (55) 


The  Euler-Lagrange  equations  for  the  adjoint  variables  \   .   are 


l  oH 


dt  ax. 


(56) 


Thus 


dAa  -  1     #    .  9fl  .  9f2    ,         1     ,   ,    ,  k     e      L%     x    U  (57) 


^T        =7<Aa-alT+Ae-a-;>+^    <  Vl+Vl+V     Sa 
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d  A   e  1  9fl  9f2  1  ?8 

"dT"     ="7   <Aa^7"     +    Ae     ^T>   +     —2    <    Vl   +    *ef2+VTI        (5&) 


dA  e  i  afi  d  f2  i  a/? 

"IT  =      "     7   (Aa-dE~     +    Ae     IT  >   +     ^2    (Aafl+Aef2  +  A  E}  TE  (59) 

B.   Control  Equation  for  Optimal  Thrust  Angle 
For  the  optimal  thrust  angle: 


9H  .0.  A  |£l  +A   5f2 


da  a  8«      e  aa  (60) 


Equations  (57),  (58),  and  (60)  are  identical  to  equations  (35),  (36),  and 
(41),  using  dE  =  —  dt.    Thus,  the  equation  for   tan  a        is  the  same 
as  equation  (50) : 


2a  AaX+  *eZ   _   N 

tan"  =   2a  X     Y  +  A  W   =  T-  <61> 

a      e 


Also,  since   t   is  not  explicit  in  equation  (52),  the  Hamiltonian 
is  a  constant: 


A  f-  +  A  f9  +  A 

H  =  — — — ~ =  CONSTANT  (62) 


C.   Transversality  Condition 

The  transversality  condition  is  now 

i 


Therefore 


[-  Hdt  +   A.dxl     +  dG  =  0  (63a) 

1   Jo 

f 

[• 


Hdt  +  A  da  +  A  de  +  AL-dE 
a      e      E 

J 


+  dt   =0  (63b) 
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Applying  the  boundary  conditions  from  equations  (38)  and  (39), 

dt  ,  da  ,  de  ,  dE  ,  da  ,  de.  -  0  (63c) 

o    o    o    o    t    r 

-  Hdtf   +  AEf   dEf  +  dtf   ■   0  (63d) 

Hence,    since     E,      is   unspecified, 

AEf  -  AE(cf>    =  °  <6Z») 

and   since      tf      is    also  unspecified, 

-  Hdtf  +  dtf  =  0  (65) 

which  yields  H  =   1  (66) 

Combining  equations    (62)  ,    (64)    and    (66) 


a   1  e   2  E 


=   1  (all    t)  (67) 


0 

and: 

£  =  Aafx   +Aef2    "  0  -  0      (t   =   tf)  (68) 

Using  equations  (67)  and  (68)  the  Euler-Lagrange  system  of  three 
equations  (57),  (58)  and  (59)  can  be  reduced  to  the  following  two  equations: 

dA,     art  /    ^£^  9f9  \ 

dt      da   \  a  d  a     e  d  a  / 

dA      ^        I  dt.  df, 

(70) 


dt     3e   I   a  ?e    Ae  de  ) 


The  solution  of  equations 

t  =  7 

3f  "  T  <72) 

f  "  f 

together  with  equations  (68),  (69)  and  (70)  is  obtained  numerically  by  choosing 
initial  values  as  follows: 

Co  =  °>  ao  '   eo  '  Eo  =  °  •  Aa(0)  '  Ae(0)  (74) 

This  requires  guessing  of  values  for  A  (0)   and  A  (0)    since  there  is 
n  a  e 

at  present  no  analytic  way  of  choosing  them.   They  must  be  chosen  so  chat  the 

numerical  integration  yields  the  required  set  of  final  values  as  follows: 

a  =  af  ,    e  =  ef,    £f  =  0  <75> 
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V.   THE  MINIMUM  E   PROBLEM 

As  mentioned  earlier,  the  minimum   E   problem  sheds  some  light 

on  the  minimum   t   problem.   We  now  will  minimize  the  final   E   (number  of 

turns  going  from  a  ,  e  ,  E   =0   to  a,.  ,  e„). 

o   o   o  f    f 


A.   Hamiltonian  and  Euler-Lagrange  Equations 


Using  equations  (30)  and  (31)  minimize   G  =  E  .   The  Hamiltonian 


is 


H  =        Aafl  +    Aef2  <76> 


The    Euler-Lagrange   equations   are 


d  A   a  9H  d  fl  9f2 

-^-  ="  TT     "     -    (    A      — r—*     +       a     — T-*  )  (77) 

dE  3a  v        a       d  a  e       9  a  v/ 


dA    e  _3H  ,  2ll  9f2. 

dE  de  (   A  a       de       +      A  e       de  (78) 

These  are  different  from  the  Euler-Lagrange  equations  for  minimum 
t   optimization. 

B.   Control  Equation  for  Optimal  Thrust  Angle 
For  the  optimum  thrust  angle: 


lii  .  0  .  »  Jll   +  x   -^  <?9> 

o<"  a   ^  «       e   d « 
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C.   Transversality  Condition 

Now,  writing  the  transversality  equation, 


1  f 


HdE  +   X  .dX. 
1   l 


+  dG  =  0 


J  o 


(80a) 


where 


dG  =  dE 


(80b) 


and  expanding  equation  (80a) 


Therefore 


o  o  o  o 

■HdE*     -   HdE,  +    \    da?   +    \    da7.  +     X   d£    + 
o  fao  af  eo 


X    d£_     +     dE,   =   0 
e      f  f 


(80c) 


-    HdE     +     dE      =  0 


(80d) 


which  yields 


Hf   "   J   =    K(l+    Xe£2    (at  E  =  E£> 


(80c) 
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VI.   TRAJECTORY  ANALYSIS 

A.   Comparison  of  Minimum   t   and  Minimum  E  Trajectories 

Equation  (80e)  is  significant;  it  means  that  any  pair  of  initial 

A's  (Aa   >Ae  )   will  automatically  generate  minimum  E   trajectories 

terminating  at  any   a   ,   e   ,   and  E   .   This  is  so  because  the  condition 

of  equation  (80e)  can  always  be  met  simply  by  scaling  the  A  's  .   Note  that 

equations  (77)  and  (78)  are  both  linear  in  the  A  's   and  that  the  optimal 

thust  angle  control  equation  (79)  relies  only  on  the  ratio  of  A  /A    and 

a   e 

not  on  their  separate  magnitudes.   Automatic  satisfaction  of  equation  (80e) 
makes  it  easy  to  calculate  minimum  E   trajectories  and  to  construct  a 
graphical  representation  of  minimum  E   transfers  on  an  e/a  plane. 

Consequently,  many  transfer  trajectories  are  generated  which  optimize 
E  ,   and  also  some  which  optimize   t  .   These  optimal  orbital  transfers  can 
be  seen  in  Figure  3,  which  is  a  graph  of  e  vs  a/a   .   It  shows  both  minimum   t 
and  minimum  E   trajectories  which  are  generated  from  the  same  initial  orbit 
around  Venus.   These  particular  transfers  concern  a  problem  incorporating 

solar-electric  orbiters  around  the  planet  Venus.   (The  thrust  acceleration 

-3    2 
used  in  the  calculations  is  5x10  '  m/s  ).   The  initial  gravitational  acceler- 

2 
ation  is  about  0.8  m/s  ;  thus  the  basic  assumption  that   A  ( (  g   is  satisfied. 

In  an  actual  case  the  value  of  thrust  acceleration  would  be  about 

-3    2  -3 

1/25  of  5x10  "  m/s  .   However,  the  value  of  5x10  "  is  used  for  numerical 

calculations  in  order  to  speed  up  the  trajectory  and  computer  times,  and  thus 

reduce  computer  costs.   Thus  all  calculated  times  must  be  scaled  up  by  a 

factor  of  25  to  achieve  a  more  realistic  mission   time. 
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Figure  4  shows  minimum  E   trajectories  only,  with  lines  of 
constant   t   and   E   superimposed.   Examination  of  the  minimum   t   and 
minimum  E   trajectories  readily  discloses  that  orbital  transfers  from  one 
point  on  the  graph  to  another  are  accomplished  in  the  least  time  by  following 
the  minimum   t  path  and  in  the  least  total  revolutions  by  following  the 
minimum  E   path.   Figure  5  shows  the  minimum   t   trajectories  superimposed 
on  the  constant   t   and  E   lines  of  Figure  4.   The  minimum  E   trajectories 
are  omitted  for  clarity's  sake. 

The  actual  calculated  orbital  elements  oscillate  over  each  turn, 
as  illustrated  by  curve  1A  in  Figure  5.   For  simplicity,  the  remaining  curves 
are  constructed  by  drawing  a  smooth  curve  through  the  values  of  the  orbital 
elements  of  each  turn  at   E  =  2n7r  (n  =  0,  1,  2,  ...).  Curve  1A  shows  the  exact 
path  of  a  minimum   t   trajectory.   It  is  obviously  an  oscillatory  type  motion, 
which  is  much  more  pronounced  here  than  it  would  normally  be  due  to  the 

increased  thrust  acceleration  (25  times  greater  than  normal).  With  the  actual 

-4    2 
thrust  of  about  2x10   m/s   the  oscillations  would  adhere  much  more  closely 

to  the  average  path  depicted  by  the  heavy  dotted  lines  -  (curve  IB)  in 

Figure  5.   The  minimum   t   trajectories,  which  are  also  minimum  fuel  trajectories 

here,  are  very  similar  to  the  minimum  fuel  trajectories  calculated  with  variable 

thrust  by  Edelbaum. 

B.   Rotation  of  Line  of  Apsides 

When  considering  optimal  orbital  transfers  between  coaxial  orbits,  one 

wishes  to  maintain  the  same  line  of  apsides.   Since  this  program  is  designed 

* 
See  Section  VIII. 
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to  change  only   a   and   e  ,   no  change  in  the  line  of  apsides  is  expected. 
A  close  watch  is  kept  of  the  change  in  apsidal  line  angle  throughout  all 
orbital  transfers  and  it  is  verified  through  calculated  results  that  such 
changes  are  small  and  relatively  insignificant.   In  some  instances  a  change  of 
up  to  one  degree  in  the  line  of  apsides  is  detected  over  a  transfer  involving 
up  to  about  50  revolutions,  i.e.,  one  degree  change  in  18,000  degrees  of 
revolution. 

Several  points  along  curves  IB  and  2  in  Figure  5  indicated  along 
with  the  flight  time  and  eccentric  anomaly  elapsed  at  each  position. 
Discussion  of  these  points  follows. 

Point  J.   Following  a  minimum  time  trajectory,  the  vehicle  arrives 
at  point  J  after  3.03  days'  time  and  37.69  radians  of  revolution.   This 
point  is  clearly  well  past  the  three-day  minimum  E   line  and  by  inter- 
polation it  is  easy  to  see  that  a  minimum  E   trajectory  would  take  about 
3.2  days  to  transfer  to  the  same  orbit.    Now,  while  the  minimum  time 
trajectory  has  gotten  to  the  same  point  in  less  time  than  the  minimum  E 
trajectory,  it  has  gone  through  more  degrees  of  revolution  to  accomplish 
this.   Thus  the  minimum  time  trajectory  has  gone  through  37.69  radians  of 
revolution;  while,  by  interpolation,  it  can  be  seen  that  the  minimum  E 
trajectory  will  reach  the  same  point  with  only  about  36.8  radians  of 
revolution.   In  summary,  a  minimum  time  transfer  takes  3.03  days  in  37.69 
radians,  while  the  minimum  E   transfer  takes  about  3.2  days  in  36.8  radians. 

Point  L.   Here  the  minimum   t   trajectory  has  reached  the  point  L 
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in  4.57  days  while  it  takes  the  minimum  E   trajectory  about  4.8  days.   How- 
ever, it  takes  the  minimum   t   trajectory  50.26  radians  of  revolution  as 
opposed  to  only  about  47  radians  of  revolution  for  the  minimum  E   trajectory, 
Point  P.   This  is  a  very  good  point  for  comparison.   The  minimum   t 
flight  time  is  4.9  days.   This  is  clearly  less  than  the  minimum  E   time, 
which  is  well  past  the  five  day  line.   The  eccentric  anomaly  elapsed  is 
46.9  radians  which  clearly  exceeds  the  value  for  the  minimum  E   trajectory 
since  point  P  is  distinctly  below  the  minimum  E   contour  line  where 
E  =  45.23  radians. 
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VII.   DISCUSSION  OF  OSCILLATORY  MOTION 

A.  Effect  of  Thrust  Acceleration  and  Lagrangian  Multipliers  on  Oscillatory  Motion 

Through  numerous  computer  runs,  it  was  verified  that  different 
thrust  accelerations  produce  essentially  the  same  mean  transfer  trajectories, 
providing  the  Lagrangian  multipliers  are  scaled  according  to  the  ratio  of 
the  accelerations  (observe  Figure  6).   A  very  low  thrust  system  produces  a 
trajectory  with  minimal  oscillatory  motion  (curve  A_) .   As  the  thrust 
acceleration  is  increased,  and  the  Lagrangian  multipliers  scaled  accordingly, 
the  oscillatory  motion  of  the  transfer  trajectory  becomes  much  more 
pronounced  (curve  A..  )  ,  while  the  mean  trajectory  remains  the  same  as  with  the 
lower  thrust  acceleration. 

It  has  also  been  verified  (Figure  7)  that  minute  changes  in  the  initial 
Lagrangian  multipliers  result  in  uniquely  different  transfer  trajectories 
which  do  not  overlap,  as  long  as  the  same  thrust  acceleration  is  used.   The 
transfer  trajectories  can  overlap  if  both  the  Lagrangian  multipliers  and  the 
thrust  acceleration  are  changed  (Figure  8) . 

Figure  9  shows  orthogonality  between  minimum  E   trajectories  and 
constant   E   path  which  must  exist  for  optimal  solutions.   This  is  a  good 
check  on  the  problem  solution.   In  this  case,  each  point  was  plotted  to  show 
^the  detailed  angular  orientation  of  the  two  curves  at  the  point  of  inter- 
section.  The  smoothed-in  path  of  trajectory   P   in  Figure  4  does  not  appear 
to  show  orthogonality,  but  close  observation,  i.e.,  blow-up  in  Figure  9,  proves 
it  to  be  so. 

B.  Appearance  of  Negative  Eccentricities 

An  interesting  and  perplexing  problem  was  encountered  while  examining 
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the  oscillatory  points  along  trajectories  at  very  small  eccentricities. 

The  problem  is  one  which  involves  explaining  the  presence  of  a  few  negative 
eccentricity  values  which  appear  during  the  calculation  of  transfers  of  low 
eccentricity.   Figures  10  and  11  illustrate  this  point.   In  the  blow-up  of 
area  I  in  Figure  11,  points  A,  B,  C  and  D  all  lie  below  the  axis  and  are 
thus  negative  eccentricities  (of  very  low  magnitude)  and  are  printed  out  as 
such  by  the  computer. 

These  negative  points  appear  to  be  in  line  with  the  normal  cycle 
of  the  oscillatory  motion.   Another  point  to  note  is  that  none  of  the 
smoothed  trajectories  ever  go  below  zero  eccentricity;  that  is  to  say,  at 
the  end  of  every  cycle  (i.e.,  for  every    2n?r     of  eccentric  anomaly)   e 
is  always  greater  than  or  equal  to  zero.   However,  to  date  there  is  no  clear 
meaning  accorded  to  these  negative  eccentricities.   Further  analysis  is 
necessary  to  determine  the  real  significance  of  this  factor. 

C.   Oscillatory  Motion  and  the  £  Function 

Also  connected  with  the  oscillatory  motion  which  accompanies  all 
spiral  trajectories  is  the  £   function  as  defined  in  equation  (40).   Through 
many  minimum   t   program  runs,  it  was  found  that  the  £  function  never  stays 
at  zero  but  always  oscillates  through  zero  once  every  revolution,  usually  at 
every  2mr    of  eccentric  anomaly.   The  smoothed  £   function  is  zero,  but  the 
minimum   t   formulation  states  that  the  £  function  must  equal  zero  at  the 
terminal  point  in  order  to  represent  an  optimal  transfer.   Therefore,  the 
transfers  are  exactly  optimal  once  every  revolution  when  the  £  function  is 
zero  and  are  near  optimal  in  between.   The  minimum  E   transfers  are  exactly 
optimal  at  all  points  along  the  trajectory  since  there  is  no  £  function 
requirement  in  that  formulation. 
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COMPARISON    OF    TRAJECTORY    OSCILLATORY   MOTION 
FOR  DIFFERENT    VALUES  OF   THRUST    ACCELERATION 


FIGURE    6 
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COMPARISON    OF   TRAJECTORY    OSCILLATORY    MOTION 
FOR    DIFFERENT   VALUES   OF  LAGRANGIAN   MULTIPLIERS 

FIGURE    7 
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BLOW-UP   OF   MINIMUM   E   TRAJECTORY   P   SHOWING 
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OCCURENCE    OF  NEGATIVE    ECCENTRICITIES 


FIGURE    10 
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VIII.   COMPUTATIONAL  COMPARISON  WITH  EDELBAUM'S  VARIABLE  THRUST  CALCULATIONS 
A.   Examining  Thrust  Angle  Schemes  for  Different  Trajectories 

An  optimal  variable  thrust  system  can  be  utilized  to  measure  the 
limit  of  best  possible  performance  with  an  optimal  fixed-thrust  system. 
This  is  so  because  of  the  extra  degree  of  freedom  added  to  the  control  system 
when  the  variable  thrust  is  incorporated,  making  it  a  more  efficient  system 
than  the  fixed-thrust  one.   Thus,  any  optimal  performance  by  the  fixed-thrust 
system  can,  at  best,  be  only  equal  to  the  optimal  variable  thrust  system 
and,  more  often  than  not,  it  will  be  less  than  the  variable  thrust  system's 
performance. 

The  following  equations  were  derived  by  Edelbaum  (References  1,  2,  3) 
for  optimal  variable  thrust  trajectories  using  an  averaging  technique  and  are 
used  here  for  comparison  with  the  fixed-thrust  performance: 


2  -.2   2    k 

V   =  2J  t   =  A   t,   =  — 
c        f       fa 


1  + 


f-*Vf  -(Vfi.i--1 


e  -  sin 


(81) 


(AV  )2  =-   A2  T2  -  ± 
c  a 


o  L 


^aV  +  2   (Ae)2 
2a/    5   1-e2 


(82) 


A  =  the  averaged  squared  variable  thrust  acceleration  over  a 
turn  ■  a  constant 


■j**--? 


(83) 


(84) 
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J_     _1_ 

+  4 

m  _       m 
f          o 

P 

P  =  POWER  =  j  |  m  j  v. 


(85a) 


(85b) 


Equation  (81)  applies  to  large  changes  in  (a,   e) ;  i.e.,  more 

than  one  revolution.   Equation  (82)  is  a  special  case  of  (81)  and  applies 

only  to  changes  within  one  turn. 

To  compare  the  two  propulsion  systems,   V   and  aV    for  fixed 

c         c 

thrust  are  both  defined  as: 

V   =  AV   =  A  x  t  (At)  (86) 

c      c  v 

As   can  be   seen   in  Table    1,    the    fixed    thrust     V        and     aV        are 

c         c 

always  larger  than  the  corresponding  variable  thrust  V   or  a  V   where: 

c        c 

-  (AV)  fixed 
1""  (AV)  variable  ^°/; 

Note  that   C    is  always  greater  than  1.   This  indicates  a  higher  cost  for 

the  fixed  thrust  system  as  expected. 

As  might  be  anticipated  in  cases  involving  near-circular  orbits 

the  thrust  magnitude  of  the  variable  thruster  is  almost  constant,  yielding 

the  same  case  as  the  fixed  thrust  system.   Here  the  value  for   C.   is  very 

close  to  one.   When  transfers  are  conducted  with  more  elliptical  orbits, 

the  variable  thrust  system   thrust  magnitude  scheme  no  longer  resembles 

the  fixed  thrust  case,  and  the  cost  is  then  less,  making  C   take  on  larger 

values . 
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B.   Maximum   A  a   and   ie   Steering  Programs 

It  is  interesting  and  instructive  to  see  how  the  thrust  a.igle  varies 
throughout  each  transfer.   Though  only  the  angle  varies  (and  not  the  magnitude 
of  the  thrust)  in  fixed  thrust  transfers  the  angle  program  agrees  closely  to 
the  one  calculated  by  Edelbaum  (See  Figure  12).   The  thrust  at  points  A  through 
E  coincide  in  both  direction  and  magnitude. 

Examining  single  turn  characteristics,  if  e   is  changing  more 
rapidly  than  a  ,   then  the  second  term  on  the  right-hand  side  of  equation  (82) 
will  dominate  the  first  term.   Similarly,  if   a    is  changing  more  rapidly  than 
e  ,   then  the  first  term  dominates  the  second  term.   Figure  12  shows  the 
thrust  schemes  for  maximum  A  a/a   and  maximum  Ae  programs.   The  last  two 
columns  in  Table  1  show  that  when  the  predominant  change  is  in   a  ,   then  the 
thrust  scheme  is  an   "a"   thrust  scheme  and  vice  versa. 

The  fact  that  Table  1  indicates  close  but  slightly  inferior 
performance  for  the  fixed  thrust  when  compared  with  the  variable  thrust  lends 
considerable  support  to  the  validity  of  both  sets  of  calculations. 
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IX.   DETAILED  ANALYSIS  OF  THE  THRUST  ANGLE  SCHEME  FOR  A  MINIMUM  E 
TRAJECTORY 

The  minimum  E   trajectory  L   from  Figure  4  was  selected  for 
examination  of  its  thrust  angle  scheme.   The  optimal,  computed  thrust  angle 
was  recorded  and  plotted  versus  eccentric  anomaly  (see  Figure  13).   The 
trust  angle  positions  were  recorded  every  90  degrees  of  eccentric  anomaly 
for  every  third  revolution. 

Figure  4  shows  that  for  the  first  2.6  days  the  thrust  program 
changed  primarily  the  semi-major  axis.   In  this  time  period   a  doubled  from 
2  to  4x10  meters  while  the  eccentricity  increased  very  slowly  from  0  to 
about  0.03.   This  part  of  the  transfer  is  depicted  in  Figure  14  where  the 
original  orbit  is   0   ,   and  after  2.6  days  the  trajectory  has  reached 
orbit   0   .   Thus  the  program  starts  out  by  increasing  the  semi-major  axis 
with  little  change  in  eccentricity;  it  uses  essentially  a  maximum  a  a   thrust 
angle  scheme  with  mostly  tangential  thrust.   This  can  be  seen  by  observing 
the  thrust  angle  scheme  for  revolutions  one  and  four  in  Figure  13. 

By  observing  trajectory   L   in  Figure  4  one  can  see  that  between 
t  =  2.6  and   t  =  about  4.6  (days)  the  thrust  program  changes  to  a  maximum 
Ae   trajectory,  while  remaining  fairly  static  in  its  semi-major  axis. 
Figure  15  depicts  this  leg  of  the  transfer  from  orbit   0   to  orbit   0   . 
These  orbits  were  chosen  at  appropriate  points  to  illustrate  the  computed 
trajectory.   A  drastic  change  in  the  thrust  angle  scheme  between  revolutions 
four  and  seven  can  be  seen  in  Figure  13. 

Figure  16  shows  the  transfer  between   t  =  4.6  and   t  =  about  11  days 
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The  semi-major  axis   a  decreases  steadily  while   e   increases.   The  engine 
is  cut  off  when  the  trajectory  reaches  orbit   0   .   At  this  point   e  =  0.998 
and   a  =  1.24x10  meters.   It  is  interesting  to  note  that  between   t  =  4.6 
and  11  days,  as  the  program  begins  to  change   a  more  rapidly,  the  thrust 
angle  pattern  (between  revolutions  10  and  27)  begins  to  revert  back  again 
toward  a  tangential  type  thrust  scheme  (Figure  13) . 
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STARTING    FROM    ORBIT   0,  THE    SPIRALING    PROCESS    HAS 
REACHED    ORBIT  0,    ATt  =  2.6    OAYS  . 


FIRST    LEG    OF    MINIMUM    E   TRAJECTORY  L 
SHOWING    MAXIMUM     Aq    STEERING     SCHEME 


FIGURE    14 
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AT  t  -  4.6   DAYS    THE   ORBIT    SPIRALING     PROCESS  HAS 
CARRIED   THE    VEHICLE      TO    03- 


SECOND    LEG     OF    MINIMUM    E    TRAYECTORY    L 
SHOWING    MAXIMUM    Ae    STEERING    SCHEME 


FIGURE    15 
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AT    1=11  DAYS,  THE    VEHICLE    HAS    BEEN 
TRANSFERRED   TO  ORBIT    04  AT  THIS   POINT  THE   ENGINE   IS 
CUT   OFF  AND   THE  VEHICLE    REMAINS  IN  THE  SPECIFIED    ORBIT    04 


FINAL   STAGE  OF    MINIMUM  E    TRAJECTORY  L 
SHOWING  ATTAINMENT  OF    DESIRED   FINAL    ORBIT    (04  ) 


FIGURElb 
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X.   MARS  ORBITAL  TRANSFER 
A.   Description  of  Problem 

In  order  to  gain  a  practical  comparison  with  single  impulse 
chemical  retro  transfers,  realistic  values  were  assigned  to  solar  electric  and 
chemical  rocket  engines  to  be  used  in  an  orbital  transfer  about  the  planet 
Mars.   An  initial  elliptical  orbit  was  chosen  with  large  eccentricity  and 
large  semi-major  axis  (when  compared  to  the  radius  of  Mars) .   The  problem 
involved  transferring  from  this  initial  orbit  into  a  much  smaller  circular 
orbit  with  radius  of  1.1  that  of  the  radius  of  Mars.   The  periapse  of  the 
initial  elliptical  orbit  is  equal  to  the  radius  of  the  desired  final  circular 
orbit. 

The  low-thrust  spiral  transfer  is  accomplished  by  turning  on  the 
solar  electric  engine  at  the  initial  ellipse  periapse,  causing  the  space 
vehicle  to  slowly  spiral  into  the  desired  final  circular  orbit  using  the 
optimized  thrust  angle  program  for  fixed  thrust.   The  retro  transfer  is 
accomplished  by  using  a  single  chemical  rocket  retro  at  the  ellipse  periapse 
causing  the  space  vehicle  to  transfer  directly  into  the  circular  orbit.   The 
amount  of  fuel  required  for  each  transfer  is  then  compared.   The  purpose  of  the 
spiral  transfer  is  to  deliver  a  larger  payload  to  the  final  orbit.   One  must 
then  decide  if  the  increase  in  payload  delivered  is  worth  the  extra  time  required 
to  make  the  transfer.   Figure  17  shows  the  two  methods  of  transfer  being 
compared,  the  upper  drawing  being  the  spiral  transfer  and  the  lower  one  depicts 
the  single  chemical  retro  impulse  transfer  at  the  ellipse  periapse. 
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B.   Initial  Data  Used 


The  values  used  in  this  comparison  were  as  follows  (Reference  8) 


Solar  Electric  Engine 


F  =  0.59756046  kfi  '  m 


v.  =  39,464.9  m/s 


Chemical  Rocket  Engine 


v.  =  2890.49  m/s 
J 


Both  Engines 


M  nn  =  initial  mass  in  ellipse  =  2049.5093  kg 
ell 


Mars 


radius  of  Mars:   3400  km 
13  3,2 


=  4.293x10   m/s 


a  ,,  -  9.734764xl07  m 
ell 

e  =  1  -  — —  =  0.95773909 

a  1 1 
ell 


C.   Chemical  Rocket  Transfer 

The  Mathematical  Development  follows: 


ve  -  velocity  at  ellipse  periapse 


-A 


•  2k    k  ' 


(88) 


=  4.6389xl03  m/s 


v   =  velocity  at  final  circular  orbit 
c 


-VF-3. 


390  m/s 


AV   =  v    -  v  =  1,248.9  m/s  (89) 

c    ep    c 
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mc 


-  mass  in  final  circular  orbit  after 
chemical  retro 


mClAn  ,,  =  e"Av/vj  =  e"1248 . 9/2890.49 
4   p  J    e 


l'  ell    c     J  -  e  (90) 

-.431 
~  e      =  .652 

So   mCi  =  (.652)(2049.5093)  =  1,340  kgs 

Amc   =  2049.5093 
-1340.0000 
709.5093  kgs 

Thus,  the  direct  chemical  rocket  insertion  requires  709.5093  kg 

of  fuel  with  an  initial  total  vehicle  mass  of  2049.5093  kg.   This  means 

that  a  total  mass  of  1,340  kg  is  delivered  into  final  orbit  by  this  method. 

D.   Optimized  Fixed-Thrust  Acceleration  Spiral  Transfer 

Next,  the  same  transfer  is  accomplished  using  the  solar  electric 
spiral  technique.   Assuming  the  same  initial  ellipse  and  total  vehicle 
mass,  the  initial  thrust  acceleration  is  determined: 

F  -4    2 

A  =  — ^—  =  2.92x10   m/s 

ell 

Although  assumed  a  constant,   A  will  actually  increase  as  the 

-4    2 
vehicle  mass  decreases  due  to  loss  of  fuel.  A  mean  value  of  about  3x10   m/s 

was  used  in  the  programs,  which  implies  a  total  anticipated  change  of  about 

6  percent  in  thrust  acceleration.   The  time  to  spiral  optimally  to  e   ,   a 

was  computed  to  be  128.685  days.   To  compute  the  fuel  expended  while  spiral ing 

into  the  final  circular  orbit  (Reference  9): 


F   =  |  mf  v.  (91) 

|*|-     f  " 

11         v . 


2 
F_        .59756046      kg. m/s 

39,464.9  m/s 


m   =    1.515xl0"5    kg/s 
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To  find  the  fuel  expended  while  spiraling,  Am   ,   multiple  |m| 
by  the  total  time  t  elapsed  during  the  spiraling  process: 


Amg  -  |m  |  t  (92) 

Now, 

t  =  128.685  =  1.11716xl07  sec 

Amg  =  (1.515xl0~5  kg/s)(1.11716xl07  s) 

=  169  kg 
E.   Results  and  Comparison  of  Data 

Thus,  the  solar  electric  spiraling  process  has  accomplished  the 
same  transfer  while  using  only  169  kg  of  fuel.   This  yields 

Initial  mass       2049.5093 

Fuel  -   169.0000 

m„  1880.5093  kg 

c2 

1880.5093  kg  of  mass  is  delivered  into  final  circular  orbit  by  the 
solar  electric  method.   Since  709.5093  kg  of  fuel  were  required  by  the  chemical 
retro  transfer,  and  only  169  kg  were  expended  during  the  solar  electric  spiral, 
the  latter  method  yields  a  savings  of  540.5093  kg  or  about  1150  pounds. 

These  results  mean  that  the  desired  final  circular  orbit  can  be 
obtained  almost  instantaneously  by  chemical  insertion,  but  with  a  cost  of 
709.5093  kg  of  fuel;  or,  the  same  orbital  transfer  can  be  accomplished  over  a 
period  of  128.685  days  with  a  cost  of  only  169  kg  of  fuel.   Thus,  if  one  chose 
the  slower  spiral  technique,  the  amount  of  payload  delivered  to  final  orbit 
could  be  increased  by  about  1150  pounds.   The  advantages  and  disadvantages  of 

*~~ 
Note  that  for  constant  thrust  magnitude,  the  total  increase  in  thrust 

acceleration  due  to  fuel  depletion  =  9  percent. 
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each  method  must  be  weighed  along  with  the  mission  objectives  to  arrive  at 
a  decision  over  which  method  should  be  used.   It  would  seem  that  the  spiral 
technique  as  described  would  be  far  better  suited  to  an  unmanned  observational 
missions  than  the  chemical  method  (Reference  10) .   It  would  seem  impractical 
to  use  the  slow  spiral  technique  on  any  manned  flights,  regardless  of  the  fuel 
saved,  but  this  does  not  preclude  the  use  of  solar  electric  engines  enroute 
to  various  destinations  as  a  thruster  boost. 
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XI.   CONCLUSIONS  AND  RECOMMENDATIONS 

The  numerical  solution  of  the  optimum  fixed-thrust  acceleration 
elliptic  spiral  technique  has  been  used  to  generate  many  different  planeto- 
centric  transfer  trajectories.   Figure  2  represents  an  infinite  number  of 
minimum  time  and  minimum  eccentric  anomaly  transfers.   It  should  be  noted  that 
any  transfer  generated  can  also  be  reversed  by  appropriate  manipulation  of 
the  Lagrangian  multipliers.   Optimized  transfers  can  be  conducted  from  any 
point  along  a  generated  trajectory  to  another  point  on  that  trajectory.   This 
is  done  by  starting  the  program  with  the  Lagrangian  multipliers  already 
printed  out  at  the  desired  point  along  the  trajectory. 

At  present  there  is  no   analytic  method  for  choosing  the  Lagrangian 
multipliers  for  particular  transfers.   By  generating  many  transfer  trajectories, 
one  can  accumulate  many  cases,  but  the  initial  process  involved  in  choosing 
the  Lagrangian  multipliers  remains  guesswork  and  experimentation.   Possibly 
by  examining  more  closely  the  transfers  produced  by  the  numerical  solutions 
used  herein  one  might  draw  some  clues  for  the  analytic  solution.   The  use  of 
reiteration  techniques  to  converge  initial  guesses  to  final  correct  values 
was  not  explored  in  this  thesis,  but  is  a  well-established  art  (Reference  11). 

As  mentioned  in  an  earlier  section,  all  transfers  were  designed  to 
be  coplanar,  and  co-axial.   As  expected,  virtually  no  rotation  of  the  line  of 
apsides  was  observed.   A  problem  noted  and  recommended  for  further  study  was 
the  occasional  production  of  very  small  negative  eccentricities  along  certain 
transfer  trajectories.   It  only  occurs  in  transfers  involving  circular  or 
near-circular  orbits  and  then  only  in  transfers  displaying  maximum  a«  thrust 
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angle  schemes. 

The  hypothetical  Mars  orbital  transfer  conducted  under  realistic 
conditions  yields  some  interesting  facts  for  consideration.   The  optimized 
spiral  technique  proves  to  be  able  to  deliver  up  to  1150  pounds  more  into  the 
desired  final  orbit  than  the  single  impulse  chemical  retro  method,  based  upon 
a  decrease  in  the  amount  of  fuel  required  to  accomplish  the  transfer.   However, 
the  spiral  transfer  requires  128  days  to  make  the  transfer,  while  the  chemical 
method  works  almost  instantaneously. 

Another  possibility  would  be  to  assume   that  no  more  payload  is 
required  on  this  particular  flight,  in  which  case  the  extra  1150  pounds 
could  be  trimmed  from  the  initial  mission  weight.   This  would  greatly  reduce 
the  fuel  required  to  launch  the  vehicle  from  earth  orbit;  or,  one  could  use 
the  same  amount  of  booster  fuel  in  which  case  the  vehicle  would  reach  the 
target  planet  much  faster  than  before  and  would  thereby  decrease  the  overall 
flight  time.  One  must  then  look  at  all  of  the  factors  involved  in  particular 
missions  to  determine  which  method  or  combination  of  methods  might  be  most 
useful. 

While  this  thesis  has  addressed  itself  only  to  the  use  of  optimized 
fixed-thrust  orbital  transfers,  it  is  quite  conceivable,  and  possibly  even 
more  advantageous  to  consider  using  the  electrostatic  engine  in  heliocentric 
space  enroute  to  the  target  planet  as  well  as  after  planetocentr ic  space  has 
been  entered.   Any  judgment  on  how  to  interface  a  heliocentric  interplanetary 
trajectory  with  planetary  orbit  insertion  will  have  to  be  based  upon  further 
study  which  must  involve  optimizing  the  entire  flight  from  earth  orbit  to  the 
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desired  final  orbit  about  the  target  body.   Figure  18  illustrates  this  type 
of  problem.  This  involves  the  incoming  parabolic  trajectory,  chemical  retro- 
insertion  into  a  very  large  eccentric  ellipse,  followed  by  the  optimized 
spiraling  transfer  technique  until  the  desired  final  orbit  is  obtained;  in 
this  case  the  final  orbit  is  a  small  circular  one. 

It  seems  obvious  that  the  slow  spiral  method  of  transferral  would 
not  be  desirable   for  manned  flights  because  of  their  long  duration,  but 
there  are  some  instances  where  such  an  alternative  might  be  preferred.   For 
unmanned  exploratory  flights,  where  time  may  not  be  the  most  significant 
factor,  this  method  offers  an  excellent  chance  to  study  the  upper  atmosphere 
of  a  given  planet  and  to  make  planetary  observations  in  a  low,  systematic 
manner.   Starting  in  a  large  elliptical  orbit,  a  spiral  trajectory  which 
steadily  decreases  in  altitude  until  a  low  circular  orbit  is  achieved  might 
be  desirable  to  record  radiation,  magnetic,  and  electrical  fields  at 
systematically  decreasing  altitudes,  along  with  other  upper  level  atmospheric 
conditions,  planetary  observations,  etc.   Orbital  changes  conducted  by  large 
impulse  chemical  retros  make  it  impossible  to  record  such  long-time  data  as 
the  transfers  occur  so  quickly.   There  are  many  factors  to  consider  when 
weighing  the  merits  and  drawbacks  of  optimized  electrostatic  propulsion 
(Reference  10) . 

The  surface  has  only  been  scratched  in  this  field.   Improvements 
lie  in  many  directions.   For  example,  Handelsman  (Reference  12)  has  put  forth 
a  formulation  for  the  introduction  of  a  switching  function  into  the  problem 
which  would  optimize  transfer  trajectories  for  minimum  fuel  and  fixed-time 
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restraints.   This  would  allow  part  of  the  trajectory  to  consist  of  coasting 
arcs  in  which  the  engine  would  be  shut  off  by  the  switching  function  in  order 
to  comply  with  the  fixed-time  restraint. 

Based  upon  the  work  done  to  date  the  following  recommendations  for 
further  study  are  made. 

1.  Investigation  of  the  occasional  occurrence  of  negative 
eccentricities  in  certain  transfer  trajectories. 

2.  Use  of  a  thrust  acceleration  which,  for  fixed-thrust  magnitude, 
increases  as  the  vehicle  mass  decreases  due  to  fuel  expenditure. 

3.  Optimization  of  the  entire  flight  trajectory  including  inter- 
facing with  the  rendezvous  and  final  orbit  insertion  in  order  to  find  the  best 
place  to  use  electrostatic  thrusting  in  the  flight  plan. 

4.  Study  the  solutions  generated  with  a  switching  function  for 
fixed -time  problems  and  compare  with  solutions  of  the  type  covered  in  this 
thesis. 

5.  Addition  of  other  orbital  changes,  such  as  rotation  of  apsidal 
line,  non-planar  transfers,  etc. 

6.  Investigation  of  an  analytic  solution  for  a  method  of  selecting 
the  correct  initial  adjoint  variables  required  for  particular  transfers 
desired.   This  would  eliminate  the  guesswork  now  required  at  the  outset  of 
every  problem.   If  this  is  not  successful,  iteration  techniques  may  be 
necessary. 
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APPENDIX  A 
Description  of  Computer  Program 

This  is  a  description  of  the  computer  program  used  for  the  numerical 
solution  of  optimal  (minimum  time  and  minimum  E)  planetocentric  spiral 
trajectories  using  fixed  power,  fixed  thrust  propulsion,  where  the  thrust 
acceleration  is  much  less  than  the  local  gravitational  acceleration.   The 
program  is  written  in  Fortran  IV  computer  language  and  incorporates  a 
straightforward  Runge-Kutta  type  integration  process.   It  was  designed  for 
use  on  the  360/30  computer  system,  but  can  easily  be  adapted  for  use  on 
other  systems.   Double  precision  was  used  in  all  program  calculations 
(Reference  13) . 

In  this  appendix  the  important  program  input  and  output  parameters 
are  described,  followed  by  the  program  itself  and  the  program  flow  chart. 

Program  Input  Parameters 

In  order  to  produce  a  trajectory  for  optimal  orbital  transfer, 
certain  initial  quantities  must  be  fed  into  the  computer.   These  can  be 
divided  basically  into  three  sets  of  values;  the  initial  orbital  conditions, 
the  two  Lagrangian  variables  (A  ,A  )>  and  the  thruster  conditions. 

The  initial  orbital  conditions  include: 

a   (initial  semi-major  axis) 
o 

e   (initial  eccentricity) 
o 

k   (mass  constant) 

t   =0  (transfer  time) 
o 

E   =0  (eccentric  anomaly) 
o 

w  =  0  (angular  position  of  line  of  apsides) 
o 
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The  thruster  conditions  include: 

A  (thrust  acceleration) 

Ml  (number  of  revolutions  desired) 

M2  (number  of  steps  per  revolution) 

The  initial  Lagrangian  variables  are: 

X  (0)gl 
a      a0 

e      eo 

The  accuracy  is  determined  by  M2  (the  number  of  steps  per  revolution 
desired).   In  other  words,  if  M2  =  50,  then  the  program  will  compute  and 
print  out  the  spacecraft's  coordinates  50  times  during  each  revolution,  giving 
the  exact  description  of  the  elliptical  orbit  which  would  be  obtained  were 
the  engine  to  be  cut  off  at  any  one  of  those  points.   However,  we  have  found 
through  experimentation  that  a  value  of  M2  ■  10  is  sufficiently  accurate  and 
generates  trajectories  more  quickly. 

For  every  initial  set  of  Lagrangian  variables,  a  different  transfer 
trajectory  will  be  generated.   The  type  of  transfers  generated  are  indicated 
in  Figure  4. 

Computer  Output  Parameters 

At  each  point  in  the  transfer  trajectory  a  number  of  instantaneous 
quantities  are  printed  out.   These  include  the  total  step  number,  which,  if 
divided  by  the  number  of  steps  per  revolution,  yields  the  number  of  revolutions 
to  that  point  (a  total  step  number  of  156  would  mean  15.6  revolutions  at 
ten  steps  per  revolution);  the  eccentric  anomaly  in  radians;  the  thrust  angle 
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in  degrees;  the  £  function,  which  must  go  through  zero  for  an  optimal  transfer 
to  take  place;  the  eccentricity;  the  change  in  eccentricity  per  radian  of 
eccentric  anomaly;  the  semi-major  axis;  the  change  in  semi-major  axis  per 
radian  of  eccentric  anomaly;  the  Lagrangian  variable  A    and  its  change  per 
radian;  the  Lagrangian  variable  A    and  its  change  per  radian;  the  time 
elapsed  since  thrust  initiation;  the  number  of  seconds  of  time  per  radian  of 
eccentric  anaomaly;  angle  of  the  line  of  apsides;  the  change  in  the  line  of 
apsides  per  radian  of  revolution. 
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Program  Listing 
CdXDQ)  MAIN 

DOUBLE    PRECISION    H ,X, Y, Z, A,GMO ,N, D, H,SIGB 

COMMON  W,X,Y,Z,A,GMU,N,D,H,SIGN 

DOUBLE  PRECISION  HI  (4)  ,  H2(4)  ,H3(4)  ,0(7)  ,0Q(7)  ,DQ(7)  , DELQ (7) , D1 , D2, 
1  T,  TT,  PR(7) 

NAMELIS7  /INPUT/  M 1 , M 2, A, GMU, Q, SIGN 

DATA  Q  /  7*0. DO  / 

SIGN  =  1.D0 

DO  1 0  M  =  1  ,  4 

H2(M)  =  1.D0 

H  1  ( M )  =  0  .  5  D  0 
10  H3(M)  =  0.5 DO 

H1  (3)  =  1.D0 

ll 3  (4)  =  1.D0 

H3  (1)  =  0.D0 

H2  (1)  =  H2(1)  /  6.  DO 

H2  (4)  =  H2  (1) 

H2  (2)  =  H2(2)  /  3.  DO 

H2(3)  =  H2(2) 
20  READ  (1, INPUT, END=200) 

IF(  Ml  )  25,  200,  25 
25  WRITE(3,40)  M1  ,  H2 ,    A,  GMU  ,  SIGN 

40  FORMAT(  1H1,  15,  •  REVOLUTIONS',  15, •  STEPS  PER  REV.   A  =  ', 
1  1PD15.8,  '  MU  =  ',D15.8  ,»   SIGN  =»,D15.8  ) 

WRITE  (3, INPUT) 

T  =  O.DO 

KMAX  =  IABS(  M1  *  M2  ) 

0(7)  =  Q{1) 

X  =  M2 

DT  =  6.2831853071795865  /  X 

KOUNT  =  0 

100  TT  =  T 

DO  101  M  =  1,  6 
PR  (M  +  1)  =  Q(M) 
QQ  (M)  =  Q(M) 

101  DELQ(M)  =  O.DO 
DELQ  (7)  =  O.DO 
QQ(7)  =  Q(7) 

PR(1)  =  T  *  57.2957795 

PR  (6)  =  PR  (6)  /  86400. DO 

PR  (7)  =  O.DO 

IF(Q(6)  .NE.0.D0.OR.Q  (7)  .NE.O.DO)  PR  (7)  =  DATAN2  (Q  (6)  ,Q  (7  )    ) 

PR  (7)  =  PR  (7)  *  57.2957795 

WRITE(3,50)  KOUNT,  T,  PR 

50  FORMAT(  1  HO,  16, 1P8D14.6  ) 
CALL  DLRIV (  T,  Q,  DQ   ) 

PR(1)  =  DATAN2(  N,D)  *  57.2957795 

PR  (2)  =  DO(1)*Q(3)  ♦  DQ(2)*Q(4)  -  DO(5) 

WRITE(3,51)  PR(1),  PR(2),   (DQ  (M)  ,  M=  1  ,  6) 

51  FORMAT  (  1H  ,    8X,  1P8D14.6  ) 
KOUNT  =  KOUNT  ♦  1 

K  =  1 

GO  TO  10  3 

102  K  =  1  ♦  K 

T  =  TT  ♦  DT  *  H3  (K) 


B-2 


103 


104 

105 
106 

107 

120 


200 


CALL  DE 
D1  =  H1 
D2  =  H2 
DO  104 
DELQ  (M) 
Q(K)  = 
IF  (  K  - 
DO  1 06 
0(H)  = 
IF (  KOU 
CALL  DE 
DO  120 
PR  (M  +  1) 

PR(1)  = 
PR  (6)  = 
PR(7)  = 
IF  (0(6) 
PR  (7)  - 
WRIT?  (3 
PR(1)  = 
PR(2)  = 
WRITS(3 
SO  TO  2 
RETURN 
END 


PIV  ( 

(K) 

(K) 

H  = 

=  D 
OQ(M 

*  ) 
M  = 
00  (K 

NT  - 
RIV  ( 

M  = 

=  0 
T  * 

PR  ( 

O.D 
.HE. 

PR  ( 
,50) 

DAT 

D0( 
,51) 
0 


T,  0,  DQ  ) 

*  07 

*  DT 
1,  V 


ELO  ( 

)  ♦ 
102 

1.  7 


M)  ♦  D2 
DO  (PI)  * 
,  105. 


*  DQ(M) 

D1 
102 


)  ♦ 

KMA 

T, 
1,  6 
(M) 

57. 

6)  / 

0 

O.DO 

7)  * 
KOU 

AN2( 

1)  *Q 

PR( 


DELO(M) 
X  )   100, 
Q,  DO  ) 


107,  100 


2  957795 
86400. DO 


.OR.Q  (7)  .NE.O. DO) 

57.2957795 
NT,  T 

N,D) 

(3)  ♦  DQ(2)*0(4)    DQ(5) 
1)  ,  PR(2)  ,  (DQ(M)  ,M=1,6) 


PR  (7)  =  DATAN2  (Q  (6)  ,  Q  (7) 


,  PR 
*  57. 2957795 
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Co)a)d  DERIV 

IMPLICIT    PEAL*8     (     A-H,     0- Z    ) 
SUBROUTINE  .DERIV     (     E,    Or     DQ    ) 
DOUBLE    PRECISION     W  ,  X,  Y  ,  Z,  A,  GM  U  ,  N  ,  D  ,  H,  S  IGN 
COMMON  W,X,Y,Z,A,GHO,H,D,H,SIGN 

DOUBLE    PRECISION     E, Q  (4)  , DO  (4)  , DSI N , DCOS , DSQRT , DABS ,S IN E,COSE, 
1    C0S2E,     T1,    T2r    T3 
SQHTAU    =    DSQRT (    0(2)     /    GMU    ) 
SINE    =    DSIN  (     E    ) 
COSE    =    DCOS (     E    ) 
C0S2E    =    COSE    *    COSE 
Z    =     1.D0    -    Q  (1)     *    Q  (1) 
Y    =    DSQRT (    DABS(    Z    )     ) 
Z    =    Z    *    SINE 
X    =    0(1)     *    SINE 

W    =    Y    *     (2. DO    *    COSE    -    Q(1)     *     (C0S2E+    1.D0)) 
N    =    2. DO    *    Q(2)     *    Q(4)     *    X    ♦    Q(3)     *    Z 
D    =    2. DO    *    Q(2)     *    Q(4)     *    Y    ♦•    Q(3)     *    W 
H    =     DSQRT  (    N*N    ♦    D*D    )     *    SIGN 
T1     =    X*N     +    Y*D 
T2    =    Z*N     +    W*D 
T3    =    A    *    Q(2)     /    GMU    /    H 
DQ(2)     =    2. DO    *    T3    *    Q  (2)     *    Q  (2)     *    T1 
DQ (1)     =    T3    *    Q(2)     *    T2 

DQ(5)     =    SQRTAU    *     (1.D0    -    Q(1)      *    COSE) 

DQCO    =    -2.D0    *    T3    *     (3. DO    *    Q(2)     *    Q(4)     *    T1    ♦    Q  (3)     *    T2    ) 
1     +     1.5D0    *    DQ(5) 
DQ(5)    =    DQ(5)     *    Q(2) 
DQ(3)     =    -T3*Q(2)  *(2.D0*Q(2)*Q(U)*  (N*SINE    -    Q(1)*D/Y)     - 

1  Q  (3)  *  (2.D0*X*N    *     (Q  (1)  *W/(Y*Y)     ♦    Y*(1.D0    +    C0S2E))*D)) 

2  -    SQRTAU    *    Q  (2)     *    COSE 

DQ(6)     -    T3*Q(2)  *  (Z*D/ (Y*  Y)  *  (2.  DO-Q  (1 )  *  (COSE  +  Q  (1 )  )  )     - 

1    Y*N*  (COSE    -    Q  (1)     )     ) 
T3    =    O.DO 

IF(Q(6)  .NE.O.DO.OR.Q  (7)  .NE.0.D0)T3    =    DATAN2  (Q  (6)  , Q  (7)  ) 
T1    =    DSIN  (T3) 
T2    -    DCOS  (T3) 

DQ(7)     =    T2*DQ(1)     -    T1    *    DQ  (6) 
DQ(6)     =    T2*DQ(6)     -    T1*DQ(1) 
RETURN 
END 


/ 

01 

sk;n  -   l  .no 

NOTE   02 

BECIN   DO   LOOP 
10   M    -    1,    4 

03 


H2  (M)    -  l.DO 
HI  (M)  -  0.5DO 


04 


H3  (M)  -  0.5DO 


NO  /  END  OF  DO 
LOOP? 


YES 


06 


HI  (3)  -  l.DO 

H3  (4)  -  l.DO 

H3  (1)  -  O.DO 

H2  (1)  • 
H2  (D/6.D0 

H2  (4)  -  H2  (I) 


H2  (2)  - 
H2  (2)/3.D0 

H2  (3)  -  H2  (2) 


03.15 
20 


08 


'READ  FROM  DEV  1 

VIA  FORMAT 

INPUT 


APPENDIX  C 
Program  Flow  Chart 


->  <- 


01  .09 


WRITE  TO  DEV  3 
VIA  FORMAT 
40 
ROM  THE  LIST 


. 

15 

T    •   O.DO 

KMAX    • 
IABS(M1*M2> 

Q    (7)    -  Q    (1) 

X    -  M2 

DT   - 

6. 
2831853071795865/X 

16 

KOUNT    -   0 

03.02 

100 

17 

— 

NOTE    IB 

s   DO   LOOP 
101M    -1.6 

02.01 

19 

PR     (M    ♦    11     -    Q    (N» 
QQ    (M>     •    Q    (M) 

101 

20 

DELQ    (M>    •    O.D> 

.01 


01.20 


02 


DELQ    (7)    -   0.DO 

QQ    (71    -  Q    (7) 

PR    (1)    - 
T*57. 2957795 

PR    (6)    - 
PR    (61/86400. DO 

PR    (71    -   O.DO 


FALSE 


TRUE 


04 


PR    (7)    - 
DATAN2(Q(6)  ,    Q(7)l 


05 


C-2 


PR  (7)  - 
PR  (7)*57. 2957795 

>- 

10 

06 

i 

PR  (1)  - 

DATAN2  (N,  D) 

*57. 2957795 

PR  (2)  - 

DQ  (1)*Q  (3)  + 

DQ  (2)*Q  (4)  -   DQ  (5) 

/WRITE  TO  DEV  3   / 
/   VIA  FORMAT  50  / 
FR»I  THE  LIST/ 

NOTE  07 

11 

LIST  -  KOUNT,  T,  PR 

/   VIA  FORMAT  51/ 

notf  12 

6 

0 
1 

DERIV 
(T,  Q,  DQ) 

LIST  -PR  (1)  , 
PR  (2) ,  (DQ(M) ,  M  - 

13 

KOUNT    -    KOUNT    +    1 
K    -    1 


/   ■■;'       / 


K    -    1    +    K 
T   -  TT   +  DT*H3    (K> 


IS 


DERIV 
(T.    Q,    DQ) 


103 


Dl    -   HI    (K)*DT 
D2    -   H2    (K1*DT 


NOTE    17 


BEGIN  DO  LOOP 
104   M   •    1,    7 


II 


DELQ    (Ml    • 
DELQ    (M)    +   D2*DQ    (M) 


104 


Q    (M)    -  QQ    (M)    + 
DQ    (M)*D1 


il£*1 


ins 


Ni'U     .'.' 


106  M   -    1.    7 


03.01    — 
106 


-M 


Q    (H>    •   QQ    (M»    ♦ 
DELQ    (Ml 


02.2 


C-3 


107 


100 


03 


DERIV 
(T.Q.DQ) 


NOTE   04 


BEGIN  DO  LOOP 
120  M  -  1,  6 


120 


05 


PR(M  +   1)  -  Q(M) 


NO       /      END   OF  DO 
LOOP? 


07 


PRO)    - 
T*57. 2957795 

PR(6)    - 
PR    (6)/864O0.DO 

PR    (7)    -   O.DO 


PRO)    - 
PR(7)*57. 2957795 


WRITE    TO  DEV 
3 
VIA    FORMAT 
50 
FROM    THE    LIST 


NOTE  1  2 


LIST   -   KOUNT,    T, 
PR 


PRO)    - 
DATAN2    (N,D) 
•57.2957795 

PRO) 
DQO)*qO)      « 
DQ(2)*q(<.)      - 
DQO) 


WRITE     TO    DEV 

s 

VIA    H'RMAT 
U 
FROM    THK    LIST 


NOTE    n 


LIST  -   r*(i) 

FRO).     (I)0(M),M 

1,6) 


E=l 


C-4 


SUBROUTINE   DF.RIV    (E.Q.DQ) 
/    DERIV / 


01 


SQRTAU   - 
DSQRT    (Q(2)/GMU) 

SINE    -   DSIN    (E) 

COSE    -   DCOS    (E) 

COS2E   -  COSE*COSE 


02 


Z    -    l.DO   - 
Q(1)*(D 

Y   - 

DSQRT(DABS(Z) 
Z    -   Z*SINE 

X   -  0(1)    *SINE 


03 


Y*(2.DO*C0SE      - 
Q<1)*(C0S2E     ♦ 
l.DO)) 

N   - 

2.DO*Q(2)*q(4) 

*X  -f  Q(3)*Z 


2.DO*Q(2)*Q(<i) 
*Y  +  Q(3)*V 


H   -   DSQRT(N*N  + 
D*D)*SIGN 

Tl   -  X*N  +   Y*D 

T2    -  Z*N  +   W*D 

T3   -  A*0(2)/CMU/H 


rug 


JU. 


DO    (2)    - 
2.DO*TJ*Q(2)*q<2) 

DQ(D    " 
T3*(2)*T2 

DQ    (5) 
SQRTAU* (l.DO   - 
Q(1)*C0SE) 


Of, 


DQ(4>    -    - 

2.DO*T3(3.DO»Q 

(2>*0(<.»*T1    ♦ 

Q(3>*T2^    ♦ 

l.JDO*DQ(5) 

DQO)    - 
DQ(5>*Q(2> 


07 


DQ(3)    -    - 

T3*Q(2) 

*(2.DO*Q(2)*Q(<i) 

*(N*SINE   - 

Q(1)*D/Y   - 

Q(3)*(2.DO*X*N  ♦ 

(Q(1)*V/(Y*Y)    ♦ 

Y*(l.DO  ♦ 

C0S2E))*D))    - 

SQRTAl'*Q(2)*COSE 


08 


DQ(6)    - 

T3*Q(2) 

*(Z*D/(Y*Y) 

*(2.DO   - 

Q(l)*(COSE  ♦ 

Qd)))    - 

Y*N*(COSE    - 

Q(l))) 

T3   -   O.DO 


10 


T3    - 

DANTAN2(Q(6) ,Q(7) 


-XL- 


Tl    -  DSIN(T3) 

T2    -   DCOS(T3) 

DQ(7>    - 

T2*DQ(1)    - 

T1»DQ(6) 


_Ll_ 


DQ(6) 
T2*DQ(6* 
T1*DQ(1) 


D-l 

APPENDIX  D 
Typical  Computer  Readout 

The  following  page  represents  a  computer  read-out  of  a  particular 
transfer  trajectory  generated  by  the  minimum  time  program  previously 
described.   It  was  generated  at  a  step  size  of  ten  per  revolution,  i.e., 
step  ten  is  the  end  of  the  first  revolution.   Read  across  from  the  step 
desired.   The  meaning  of  each  column  is  indicated  at  its  head. 

When  reading  across  one  will  notice  two  lines  of  figures.   In 
each  case  the  second  figure  represents  the  change  in  value  per  radian  of 
eccentric  anomaly  as  the  trajectory  generates,  with  the  exception  that  the 
second  figure  in  the  second  column  is  the  £  function  which  must  go  through 
zero  for  optimality.   On  this  page  of  sample  read-out  the  £  function  does  not 
yet  go  through  zero.   But  later  on  in  the  trajectory  it  begins  to  oscillate 
through  zero  with  every  revolution. 


D-2 


cr- 

t>  » 

©  o 
o-  — 
<r  if 


cr> 

a- 
cr- 

»n 

r- 

o 

IT   CO 

o 

£ 

© 

O         S 

o 

f-        O 
r»<        — 

a 
t 
o 

or 

o 

r. 

o 
go 

00 
•*    00 

o 
~    1 

0  c 

1  r*> 
©   rr 

•t  o 

© 
O 

c 

o 
I 

© 
•o 

rr> 
If. 

or 
C  oc 

r-  if 

c 

?e 

c  — 
O  4 

O 

©     1 

c  o 

C  If 

41  — 

O 

c    1 

0  0 

c  c 
■*  c 

O 
©    1 
0  c 

r- 
C    «r 
0  or 

c 

C    1 

0  © 

vT 
C   r" 
4)  rr 

j 

C 

o    1 

0  c 

c 

C  o- 

4     — 

r 
C 

e 
4 

4 
O 
1 

c 
p- 

C 

1 

c 

0 
0 
1 

r 

■ 

- 

c 
0 

r 

3 

0 

1 

r 

*• 

4- 

c 
c 

c 

«c 
c 

1 

f 
'* 

-* 
4 

r 
C 

«?  C 

'      4 
C      - 

»■  ■y 
<  - 

■t 

-o 

^   f*\ 

m 

* 

CO 

D 

O  rs» 

IT    O 

or 

p-» 

0-  rs> 

r-  m 

O  rvj 

If  * 

r-  r- 

—  O 

O-    4) 

i- 

pp 

T 

•>- 

v^ 

$E 

CD 

1*1           IT 

r  i 

•4   f- 

r-  r- 

-- 

r- 

r-    mi 

r-  4 

4.  oc 

r/j  if 

^  © 

©    rv. 

a-    pr 

p^ 

p. 

.», 

c 

, 

J 

o 

>r 

cr 

o 

<x 

r-  o 

od  rsi 

cr 

pi 

\/\  on 

pi  pi 

O  r« 

OC   4 

p-  r" 

O   m* 

rr    C 

< 

c 

ot 

_ 

-^ 

*T 

' 

■ 

•c     • 

IT 

f 

IT 

• 

mi      • 

f-    • 

tr 

• 

4     • 

m     • 

pp     « 

4      • 

0>     . 

«C       . 

41      • 

4 

• 

r*. 

# 

XJ 

1 

c*     • 

4  p- 

o  o 

»r  i/v 

o- 

*. 

a 

>r 

If    CD 

■r  cc 

mi 

■c 

r-  rs 

mi  pi 

rr  r- 

cr  — 

rs  0 

If    — 

4    pi 

4 

►- 

W 

«. 

^ 

'S. 

r. 

•  i 

• 

• 

I 

•     1 

• 

• 

•     1 

•    1 

•     | 

• 

1 

| 

o 

O 

p- 

— 

* 

rs 

fM 

* 

«9 

4 

•«■ 

rs 

pi 

m 

fT 

wm 

»*^ 

4 

r*l 

*> 

«^ 

i»i 

m 

1    f 

m 

r*. 

rr 

1    ^ 

1  •" 

1     " 

|P\ 

ip- 

fp 

fT 

p 

m 

o 

© 

o 

o 

O 

c 

o 

O 

o 

O 

0 

c 

O 

0 

0 

c 

O 

c 

c 

m 

CM 

fNJ 

rg 

rs 

rM 

rs 

rs 

rs 

rv. 

rs 

rs 

pp 

pp- 

mm 

_^ 

#_- 

a 

o  o 

O  O         O 

c 

©  C 

o  C 

o 

© 

C   C 

O  © 

?8 

O  © 

0  c 

©  C 

c  c 

0 

c 

c 

c 

c 

C 

r 

c 

c  c 

o 

m 

1  ►- 

1 

1^       1 

* 

1        — 

I    -c 

j 

e> 

1    4 

1    4 

1    — 

1     re 

1    O 

1   c 

1 

4 

r* 

1 

i^- 

rs. 

1    c 

c 

cr- 

C or. 

c 

o- 

c 

c  c 

c  o 

c 

wm 

c  cr 

C    <-> 

c  c 

C  — 

C    O 

C    — 

f      4 

c 

4 

c 

«r 

f 

r 

4T 

r  4 

C 

>, 

oc. 

o-  o- 

>» 

h        m 

~r 

c^  »*r» 

f-    IT 

ir 

rr\ 

•c  o- 

IT  «- 

rs  rs 

*  — 

4.    rr 

p-   p- 

IS   c 

.f 

c 

f^ 

r . 

^ 

4 

V 

cr     4 

c 

• 

m 

rvj  rr» 

-- 

* 

a 

-r 

00  ^ 

C   * 

r- 

■*• 

i  rn 

4  CN 

rs  ". 

4  f> 

rs  r" 

r-  4 

r-   -r 

4 

If 

-c 

tr 

-C 

*r 

47 

r  - 

c 

c 

-o 

o 

o  o 

r- 

C 

•» 

© 

o  o 

r\j  o 

rr» 

O 

r-  c 

o-  c 

r-  O 

If  C 

or  C 

*  0 

4   C 

c 

C 

c 

cr 

c 

rv. 

c 

pT     C 

c 

*--■ 

pp 

on  _ 

ir 

cr 

•■ 

*•  *rt 

sr  m* 

op 

pi 

p-  pi 

rr-  pi 

41  — 

0D   — 

c  — 

m"    pi 

ip-   p- 

4 

4 

tf 

tf 



4    _ 

o 

• 

• 

<N      . 

C 

• 

• 

»*      • 

-*    • 

— 

• 

rs    • 

rs    • 

rs     • 

rs     • 

rr      • 

in     • 

*      • 

IS 

1 

c 

• 

a 

• 

W 

1 

o 

«-» 

Opp 

o  — 

* 

<r 

fM  mi 

o  — 

or 

•p. 

•c  — 

4  — 

rs  p- 

O  — 

or  — 

4)  — 

C   — 

PI 

pp 

tv 

,— 

f^ 

„,! 

pt 

_„ 

4   p" 

© 

• 

• 

• 

• 

t 

• 

• 

• 

t 

• 

o 

o 

■ 

•■■ 

r«j 

m 

■«■ 

-r 

IT 

•O 

r*- 

or 

or 

cr 

pp 

— 

.M 

^ 

_ 

_ 

© 

•f 

■r 

<r 

* 

J- 

»f 

■♦ 

•r 

•J 

-J 

4 

4 

* 

.*• 

4 

•* 

-# 

^ 

^ 

o 

1    • 

o 

o 

l\l     1 

o 

«V     1 

»Si 

o 

i 

f\. 

o 

i 

r,? 

r.C 

•M 

? 

c 
rs    | 

c 
rs   I 

O 
rs    1 

O 
rs    1 

© 
rs    1 

c 

c 
rv    1 

IS) 

C 
1 

f^ 

r 
1 

»^v 

<? 

fs_ 

' 

C 
^.     1 

pi 

• 

O  Q 

c  o 

o 

O        c 

c 

o  c 

c  c 

O 

C 

©  C 

©  C 

<?£ 

?5 

e  r 

©  r 

r  r 

e 

r- 

0 

c 

c 

c 

C 

r 

E    ' 

Q 

1    -O 

{     IT, 

1 

">       1 

■r 

1  *> 

i   tr* 

1 

»»\ 

1    O 

1   t» 

1    4 

1   4 

1    4 

r 

p^. 

1 

1-- 

CT 

1 

c 

.    p- 

■ 

c  o 

C  <n 

c 

<\ 

C 

J- 

C  .e 

C  rs> 

C    fr\ 

C   4 

C  IT 

C    vT 

c  e> 

c  e 

c  — 

©  if 

fw 

c 

^ 

c 

g 

p- 

c  c 

■ 

O  r- 

oc  o 

mi 

•» 

•c 

p*. 

—  oc 

rr-  rsj 

mi 

r- 

<r  c 

f»   tr\ 

cm  a 

r<-  rs 

<r  « 

If    IS 

pp  rr 

4 

'. 

— 

•_ 

4*-- 

Ml 

»r 

rv 

.    - 

ST 

a 

1 

•      ©  -i 

<£  (<> 

K" 

r- 

^ 

m> 

IT   41 

oc  r- 

rs. 

cr 

0-  or 

0  c 

4;  rr 

pp)  ► 

C    tf 

a-  ■* 

rr 

rr- 

T 

4- 

1*" 

4» 

c 

m- 

•"     4 

o 

• 

-< 

o  r- 

C  M 

•• 

ir 

fw 

^ 

0"  -i 

<M  r- 

•o 

pi 

r-  <. 

r~  4 

*  c 

ss 

m  O 

C   J- 

ip-   43 

PI 

rr 

t/ 

c 

c 

m 

1^- 

4 

■     p 

o  0> 

ir  © 

C 

m 

r- 

rv 

•»>  IM 

C  pi 

* 

rsi 

r*  rs 

C  — 

If.  0 

if  C 

C  — 

»-  rs 

4 

r- 

c 

I*-. 

rv 

»-• 

p 

©    pp 

tjfl 

r*   • 

c    ■ 

•0 

< 

CC 

• 

p-     • 

*    • 

~c 

• 

a    . 

pi     • 

*    • 

r»     • 

f>     • 

rs    • 

vT      • 

p» 

• 

c 

• 

r\. 

• 

*r 

• 

or     . 

•»  m. 

K  -f 

r- 

«r 

r» 

•r 

•=  ♦ 

oc  * 

a 

^> 

oc  * 

<r  * 

0  •» 

C  r<- 

c  * 

©  4 

C    ♦ 

C 

■T 

■■ 

-» 

•«- 

4- 

•— 

4 

—  ■*• 

• 

• 

• 

• 

• 

• 

t 

■ 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

• 

m 

e* 

.* 

•M 

•-* 

ml 

— 

•* 

p-< 

pi 

pi 

pi 

pp 

pi 

rs 

rs 

rs 

f\ 

f>w 

r\ 

rv. 

PI 

c 

m 

r\j 

r 

r\ 

«\i 

r- 

f*. 

rv 

rs 

rs 

r~ 

ISl 

rs 

rs 

rs 

(■»■ 

rv. 

r 

r 

« 

O 

o 

o 

O 

o 

© 

© 

C 

© 

0 

© 

O 

O 

O 

c 

c 

C 

C 

C 

c 

o 

c 

* 

ir 

ir 

IT 

%f> 

IT 

IT 

IT 

ir 

if 

if 

if 

if 

4 

«r 

*r 

IT 

j- 

4 

o 

o> 

• 

o  o 

O  Q 

r 

B 

O 

c 

o  c 

o  c 

C 

c 

c  c 

O  © 

c  c 

c  c 

©  c 

c  c 

c  © 

©  c 

C 

c 

c 

r 

c 

C 

C    C 

o 

0< 

c 

<•> 

o 

m 

tr. 

•* 

0- 

rA 

«r 

CD 

4 

r- 

c 

4) 

r- 

4 

ec 

K 

p 

4 

c 

c 

©-< 

C  u-> 

c 

in 

C 

r- 

C  (*-. 

t-  m 

r- 

pi 

r»  .» 

rv  r- 

C  rv 

c  * 

C   or 

C    4 

C    4 

© 

4 

c 

O 

C 

P* 

c 

c 

r   rs 

o 

01 

• 

©  o 

^  © 

-f 

f- 

rv. 

r 

•f  0> 

^  -- 

rsi 

O 

—  j» 

oc  ir 

or  y 

vf      -• 

©   •»■ 

If  p- 

l»    C 

4 

4 

« 

*? 

oT 

4 

H 

a 

rs    4 

© 

en 

C  r- 

O"  rn 

o 

© 

O 

<■ 

F>    — 

r-  r~ 

r- 

fw 

-»   j- 

<n  C- 

4  0 

*      ►- 

rr   If 

«»•  r- 

^    ~ 

pi 

^. 

-— 

i»- 

.SB. 

pp 

a 

4 

or  r- 

r~ 

en 

1 

»       c  ""1 

or  «M 

ri 

— 

O 

0D  C 

0*  r»- 

O 

a 

*   C 

or  rv- 

O  rvj 

CT   rr\ 

a-  1-' 

©  •* 

c  c 

a: 

c 

r 

I-" 

r 

C" 

fc 

r 

p    pp 

cr- 

en 

■« 

o  o 

m  ■*■ 

© 

IM 

c 

r  j 

w  <• 

e  c 

43 

r- 

0-  rs 

C  rs. 

4>  J- 

©  © 

rr  -r 

©  rv. 

©  fs 

rr. 

.!• 

c 

c 

<c 

IP- 

C* 

a 

C*    rv. 

fs> 

8 

c    • 

c    • 

o- 

■ 

a 

• 

0>     « 

o    • 

© 

• 

o    . 

C     • 

C     • 

cr    • 

o-    • 

a    • 

C      • 

c 

• 

r 

• 

r 

• 

r 

• 

C      • 

• 

«r  — i 

.»    CO 

-I 

f-1 

4 

(H 

■♦  oo 

■«•  -< 

IT 

rr- 

If  rr. 

»T   rr- 

if  or 

4  -< 

4   IT 

4  m 

*   !»■ 

4 

or 

<r 

•— 

»r 

c 

«r 

p- 

4    p- 

* 

rvl 

•  i 

•    I 

• 

1 

• 

• 

• 

• 

• 

•  1 

•    1 

•   1 

•    1 

•    1 

• 

• 

• 

• 

• 

•     1 

s 

tn    ' 

«-, 

r-> 

f 

l»* 

■*> 

PI 

1* 

f 

rr 

IT. 

pr. 

ir- 

m 

pr 

r» 

^ 

•*■ 

r- 

«S) 

<M 

rj 

r\ 

r\i 

rsi 

fsi 

rv 

rs 

rv 

rs 

rv. 

IS 

rv 

rs 

IV 

rv 

- 

p 

* 

4 

O 

a 

O 

o 

© 

o 

© 

o 

© 

© 

O 

© 

C 

© 

O 

c 

C 

C 

O 

■e 

<■ 

< 

-c 

41 

41 

4 

4! 

•c 

< 

4 

•c 

4 

4 

4 

V 

^, 

£ 

~ 

3 

o  a 

©  n 

©  c 

©  c 

o  © 

©  o 

o  c. 

©   C 

o  © 

©  C 

0  c. 

c  © 

©  c 

c  c 

a 

c 

r 

r 

c 

r 

t- 

P 

C     ' 

X 

« 

n 

r- 

r^ 

X 

<c 

rr. 

-r 

c 

rr. 

r-. 

p* 

rr 

c 

1*- 

— 

3 

p- 

*- 

rf 

3 

©  on 

C  M 

c 

m 

c 

ir 

C  r~ 

c  <y 

c 

r^ 

C      vf 

c  c 

C  ec 

S5 

C    V 

C   <r 

©    © 

c 

or 

c 

*— 

c 

4 

r 

' 

c     - 

-r 

r 

C  «f> 

r-  im 

IT  o 

r- 

c 

r»-  O 

r-  a 

pi 

r- 

—  4 

rs  rs 

if  — 

rr  m 

<r  4; 

rr  if 

pp 

4 

r 

•— 

T 

4 

c 

vT 

r  c 

o 

o 

O  -o 

^   -O 

rr 

r- 

J- 

oc  *+ 

«»   r>- 

p— 

<. 

<s  rv. 

r-  r- 

•c  4 

©  cr 

If    4 

*  4. 

c    •* 

rv 

u- 

C 

<— 

*r 

r 

^ 

~ 

^  - 

© 

• 

■ 

o  •#■ 

«r  >r 

or  r- 

o- 

<  --> 

rvj  sf 

oc 

rr 

IT  r- 

4(  4 

C    4 

«■  4 

©    r- 

rr   © 

4    r 

rv 

e 

»>-■ 

fcT 

r-- 

a 

c 

fc 

r      p 

fl 

O  — l 

c  ff5 

O 

o 

f* 

o  •*■ 

©  — 

c- 

4" 

o-  ir 

O  rs 

0  or 

0  — 1 

O  on 

p-  r 

pp     I* 

pp 

rr. 

C 

«■ 

*r 

r 

c 

p- 

C     l» 

© 

1  "* 

a 

>»    • 

^    • 

* 

I 

* 

■ 

■♦r    • 

-r     • 

m 

• 

rr      . 

rr      m 

vT      • 

4    • 

^      • 

4     • 

^»    • 

J 

• 

-" 

' 

o- 

• 

* 

• 

4       • 

c 

e 

• 

• 

• 

. 

1 

•  i 

•    1 

• 

1 

•   1 

• 

• 

• 

, 

, 

•    1 

. 

1 

• 

| 

• 

. 

, 

• 

4 

c 

• 

•»> 

r 

r- 

r» 

r" 

fr 

f*> 

**■ 

c- 

fp- 

r* 

rr 

rr 

p 

ft 

t* 

»-- 

" 

p 

f 

c 

© 

■«■ 

•r 

-» 

-r 

<f 

•r 

>r 

vj 

J 

»T 

•jr 

*» 

•» 

*T 

4 

sT 

4 

h 

- 

IN 

c 

o 

a 

c 

c 

O 

c 

© 

o 

© 

4? 

© 

C- 

© 

r 

e 

C 

C 

© 

© 

«r, 

en 

•> 

1 

•j-   i 

^ 

1 

* 

1 

•r  i 

*   1 

4 

i 

v»       1 

4     1 

4     1 

rr     1 

©cJ 

rr     1 

rr 

1 

** 

1 

•^ 

1 

F 

1 

- 

» 

en 

pp 

a 

c  o 

C 

r 

C 

c 

o  c 

o  c< 

c 

c 

c  c 

O  C 

c  c 

c  c 

©   C. 

C  f 

c 

c 

c 

r 

e 

t. 

C 

• 

© 

o 

-< 

1     K- 

• 

* 

1 

© 

<.£ 

i 

IT- 

1  »<-. 

1    4 

1     4 

1     4 

1     or 

1      © 

1    c 

1 

4 

1 

•» 

1 

1 

* 

•V, 

o> 

1 

«o 

c  o< 

c 

C 

C 

|T 

rsi 

C  a 

c  <•- 

r    r>- 

r  isi 

C    C 

r  c- 

r-      If 

C 

pp 

c 

^ 

c 

— 

. 

r  — 

0" 

c 

oa 

^  on 

-- 

e 

>C 

•T 

oc  O 

»»■  c 

IT 

pi 

-»  vT 

—  © 

•c  cr 

If    4) 

C    -i 

►-    4 

p-    P- 

rf 

4 

l-s. 

r 

•V 

< 

*^ 

p- 

-  - 

rv 

© 

*< 

a 

IT  — < 

C 

<r 

'       c 

«T 

#r>  en 

r-  •-• 

<7 

r— 

tr  ^ 

rr  o 

C    rj- 

CD    O 

•0  c 

a-  p-. 

a    1 

pi 

4 

**-. 

i 

•— 

W 

c 

» 

f  ' 

• 

r 

•H 

©  j- 

~r 

a- 

•       u- 

p\ 

©  f' 

or  -^ 

«r 

rr* 

r~  -i 

or  w 

rr   rr 

vt     © 

4-    y 

rs   J 

0    C 

a 

m 

— 

•—• 

iV 

1* 

-T 

4 

^ 

■» 

O'- 

c 

r-  ^1 

•M 

c 

>       * 

C 

*  o 

•C   0" 

a 

*< 

4  c 

c  © 

—  4. 

m  er- 

»    4 

—  O 

0   c 

rv 

4) 

K 

0 

— 

-T 

J 

C- 

r      C 

m 

r" 

is 

er 
0 

1 

©  — 

.»     « 

O 

■        cr 

fvj     < 

4     • 

o 

4      • 

rr      • 

a     • 

rs     • 

IS       • 

"•     * 

*     * 

V 

* 

a 

tf 

C 

3  - 

rr 

C 

• 

i 

• 

• 

• 

I 

• 

• 

• 

• 

• 

• 

1 

• 

• 

• 

• 

• 

• 

• 

© 

0 

© 

mt 

■^ 

rv. 

r" 

•r 

tr 

.r 

r~ 

tv- 

O 

— 

p- 

^ 

-p 

— . 

— 

— 

— 

> 

rs 

a 

,-^ 

f*- 

n 

«• 

1 

m- 

e 

c. 

rr\ 

r*. 

f» 

rp> 

rr 

pp 

|» 

pr 

r» 

" 

f 

" 

u* 

c 

cr 

m 

c 

o 

C 

) 

C 

C 

c 

© 

© 

O 

O 

© 

© 

C 

c 

r 

c 

C 

c 

C 

a 

c 

<v 

prf 

rs. 

fVJ 

rs 

r^ 

rv 

rv 

rs 

rs 

r 

IS 

r 

r 

1* 

IN 

fv 

r 

c 

<u 

c 

o  o 

©  c 

1         C 

r 

O  C 

C    C 

O 

cu 

c  r 

r  C 

c  c 

C  C 

©   C 

C    r 

c    c 

C 

C 

c 

c 

C 

1 

f 

( 

or 

H 

cr 

01 

<•■ 

u> 

o 

IT 

t              rs« 

? 

j 

^ 

If 

O" 

if 

vf 

p 

0 

•" 

— 

*• 

*■ 

UJ 

« 

©■ 

41 

0> 

r*  «N 

c 

•t 

3          C 

f\ 

C  r- 

r  sr 

c 

fp- 

c  c 

C  oc 

C    if 

©  on 

r   r- 

©    f 

r   j- 

C 

•*» 

r 

■^ 

r 

c 

r 

a 

I 

>  c 

-o 

•- 

C  *- 

c 

r 

4           C 

r 

c  c 

©  r~ 

c 

o- 

C    «Si 

©  ^ 

C   rp. 

c  — 

C    p 

C   or 

C     4 

<- 

r 

c 

•>■ 

4 

f 

C     — 

c 

r- 

4 

c  -»■ 

c 

a 

1       c 

C  or 

C  © 

r- 

a. 

c  — 

r  -f 

C   rvj 

C     IM 

C    r- 

©    4 

C    r 

C 

—- 

c 

«• 

c 

■~ 

c 

■ 

</> 

m 

„ 

UJ 

f\ 

c  n 

c 

u 

>       © 

K 

©   IT 

r  *-■ 

r- 

IT 

r    r- 

O  IT 

C    <*. 

t~   rs 

C    m 

r  4 

C    1- 

C 

j 

r 

0 

c 

•J 

r 

p- 

r    . 

c 

• 

c 

c  © 

© 

r 

J       c 

C 

>       c  C 

c  c 

C 

© 

C    C 

©  © 

C    © 

C    O 

C    C 

c  © 

c   c 

C 

e 

c 

c 

c 

c 

c 

l 

u. 

C 

1 

c    • 

© 

■       oc 

>       >f 

>       c    • 

4 

• 

rs    • 

cr     • 

«r     • 

c     • 

4       • 

r»      • 

a      • 

4 

• 

c 

• 

4 

• 

f 

• 

■ 

p- 

c  - 

tr. 

» 

•  1 

•  i 

i 

1 

1 

•   1 

•    1 

1 

•    i 

•    I 

•    1 

•    1 

•     1 

■    1 

1 

1 

1 

e 

r- 

^> 

m 

rs 

rs 

rv. 

pr 

rr 

rr 

J 

■# 

4 

>r 

V 

- 

4 

c 

■4 

f 

1             r\ 

i             rs 

rs< 

rvj 

p- 

—* 

4 

— 

-p 

r 

rs. 

fv 

»v 

P" 

— 

pp 

© 

r 

> 

r 

>             C 

»             © 

© 

© 

© 

© 

© 

r 

e 

© 

*r 

c 

r 

r 

r- 

c 

C 

© 

C 

c 

c 

© 

© 

C     1 

c 

c 

© 

C 

c 

— 

— 

— 

cr 

C 

. 

©  C 

c 

r 

3         C 

C 

o  c 

c  r 

O 

C 

©  r 

©  r 

©  r 

c  c 

C    f 

©  r 

c  c 

c 

c 

c 

r 

c 

C 

r 

r 

C    C 

■ 

n 

0 

I    — 

r- 

pi 

J 

© 

J 

o- 

d 

rs 

pp 

© 

* 

1- 

p-- 

r 

«/ 

. 

~ 

l/i 

«v 

i  rr 

r 

a) 

r  n 

e 

r 

C 

h 

C  a 

>       c  © 

c 

sT 

C   rs 

C  4 

c   — 

r  © 

r  © 

C     IS 

r  - 

. 

a 

i. 

r 

r 

*i* 

• 

1      — 

7 

J 

c 

r» 

c 

IT    — 

C 

V 

fT 

•»  \e 

r    O 

M 

en 

r  u" 

or  r- 

4    © 

if  f 

P»    r 

r    r 

C    " 

r 

c 

r 

of 

— 

- 

c 

►  r- 

C 

a 

a  ^j 

r- 

r 

3         <T 

n 

-,        r-  j 

P  c 

^~ 

rs 

rv  *•- 

vf  rs 

4     vf 

a-  a 

C    if 

rv    a 

J    r- 

4 

J" 

1*- 

c 

— 

*\ 

& 

n 

*     . 

c 

C 

o 

""* 

—  ^ 

4 

r 

P 

r< 

rs  4 

IT    C 

c 

j 

r\  rs 

<r   or 

a    pi 

pi  n 

■    c 

0     _ 

—  r 

V* 

r- 

1** 

c 

1* 

** 

-■ 

. 

"    f 

p_ 

e 

)  B 

V 

)          --   ,o 

■r 

- 

>r 

r 

)        c  «r 

»        —  on 

a 

«- 

or  c 

jv  rs 

4    4 

r-    C 

—  4> 

r*    rv 

0   c 

» 

** 

«-* 

0 

. 

w* 

■ 

c 

r  — 

D 

c 

e 

c 

1        oc    • 

ir 

•        or 

■        •- 

•        ,»     . 

4 

• 

C     • 

rs     • 

If       • 

V       • 

p-     • 

m     • 

4     • 

C 

• 

*n 

• 

c 

• 

- 

• 

r         • 

•  B 

c 

^c  c 

1        fv.  mi. 

^ 

< 

r        oc 

-       «r  - 

(■ 

p< 

m    p- 

C  4 

4    -* 

is  cr 

C    — 

if    ■* 

—     m- 

pp 

— 

«* 

— 

c 

•■ 

c 

• 

—     ^ 

c 

c 

0 

> 

• 

a 

UIC 

•c 

IS) 

r9 

|i 

4 

«r 

W 

4 

4 

r- 

| 

a 

.- 

■~" 

— 

— 

u 

cr 

c 

r« 

r- 

r> 
C 

c 
c 

VI  © 

— 

rs 

i* 

* 

<f 

< 

r- 

o 

O 

© 

p, 

rv. 

P 

J 

* 

4 

r- 

• 

© 

~~ 

P 

cr 

a. 

_- 

— 

— 

— 

— 

— 

— 

— 

— 

e 

e 

-J" 

c 

e 

•    en 

C 

7    . 

P" 

j 

p- 

p-    p 

c 

i 

thesB1 13 

Numerical  solutions  of  optimal  fixed-thr 


1 

111 

ill 

III  I II  ill 

j 
l|    ]  I  ||| | 

■n 

3  2768  001  91122  5 

DUDLEY  KNOX  LIBRARY             r 

